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Chapter  1.  Introduction 


1.1  Introduction 

This  thesis  reports  the  results  of  an  investigation 
of  the  damping  properties  of  porcelain  enamel 
coatings.  The  aim  of  the  investigation  was  to 
determine  the  material  parameters  which  govern  the 
damping  behavior  of  porcelain  enamel  coatings. 

Damping  is  the  term  used  in  vibration  and 
shock  analysis  to  describe  any  mechanism  that 
dissipates  energy.  Damping  in  a material  keeps 
the  amplitude  of  vibration  of  a system  of  which  the 
material  is  a part,  under  control,  thus  decreasing 
the  possibility  of  failure  due  to  fatigue  or  wear. 

Damping  measurements  have  been  used  as  an 
investigatory  tool  to  aid  in  understanding  the 
fundamental  atomic  behavior  of  materials. 

Material  damping  and  the  develooment  of 
methods  of  damping  measurement  and  analysis  have 
been  dealt  with  by  Lazan  [1] . A review  of  system 
damping  and  an  extensive  bibliography  is  given  by 
Plunkett  [2].  Plunkett  [31,  Entwistle  [4]  and 
Bradfield  [5] , review  the  methods  of  damping 
measurement  that  have  been  developed  and  their 
limitations . 

If  there  is  not  enough  inherent  damping  in  a 
material,  che  system  damping  may  be  increased  by 


coating  some  parts  with  a material  with  a high  damping 
energy.  Viscoelastic  polymers  are  often  used  at  room 
temperatures.  In  this  thesis,  the  properties  of 
porcelain  enamels  as  coating  materials  at  elevated 
temperatures  are  examined.  Porcelain  enamels  are 
hard,  inorganic,  giassy  materials,  widely  used  as 
protective  coatings  on  metal  parts  exposed  to 
corrosive  and  erosive  atmospheres.  Andrews  [6]  gives 
a description  of  porcelain  enamels,  their  composition, 
properties,  uses  and  application  techniques.  They 
soften  over  a wide  temperature  range  in  which  they 
behave  anelastically . For  some  time,  it  has  been 
known  that  enamel  coatings,  in  their  softening 
temperature  range,  may  increase  the  fatigue  life  and 
vibratory  damping  of  a coated  member,  but  no  study 
has  been  made  of  the  governing  parameters. 

A cantilever  beam  apparatus  was  designed  by 
B.  J.  Lazan  to  test  porcelain  enamels,  and  the 
preliminary  tests  which  confirmed  their  damping 
capability  at  high  temperatures  were  performed  at 
Airesearch  Manufacturing  Company  of  Arizona  [7] . 

Lee  [8]  modified  the  apparatus  and  conducted  tests 
on  porcelain  enamels  at  a frequency  of  32  Hz.  Lee's 
apparatus  has  been  further  modified  during  the 
study  reported  in  this  thesis  and  two  torsion  devices 
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were  developed  to  increase  the  frequency  range  and 
investigate  the  effects  of  stress  state;  all  three 
apparatuses  are  described  in  later  sections  of  this 
thesis.  An  important  constraint  which  controlled 
the  design  of  these  pieces  of  equipment  was  that 
the  material  under  test  had  to  be  applied  as  a thin 
coating  on  a base  metal. 

Porcelain  enamels  have  an  atomic  arrangement 
similar  to  that  of  glass.  There  is  very  little 
information  available  about  the  structure  of 
porcelain  enamels.  However,  a great  deal  of  work 
has  been  done  on  glass,  and  an  extensive  literature 
exists.  In  1932,  Zachariasen  [9]  proposed  the 
random  network  structure  of  glass  which  is  now 
commonly  accepted.  This  was  later  confirmed  by 
Warren  et  al  [10] . Morey  [11]  and  Jones  [12] 
describe  the  physical  and  chemical  properties  of 
glass.  Fitzgerald  et  al  [13]  and  Spinner  [14] 
studied  the  temperature  variation  of  the  elastic 
moduli  of  some  glasses.  McCormick  [15]  developed 
a torsion  pendulum  apparatus  to  measure  the 
viscoelastic  properties  of  glass  as  a function  of 
temperature.  Forry  [16],  Kirby  [17],  Moore  and 
Day  [18],  and  Copley  and  Oakley  [19]  have  analyzed 
the  low  temperature  internal  friction  peak  in 
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various  silicate  glasses.  However,  except  for  Kadota 
et  al  [20]  , no  work  reports  the  damping  or  modulus 
properties  of  glass  in  the  softening  temperature 
range.  Most  investigations  have  been  concerned  with 
the  damping  due  to  ionic  diffusion  at  lower  temper- 
atures. This  damping  is  very  small  and  cannot  be 
used  for  engineering  purposes.  An  Arrhenius 
relation  [21]  between  time  and  temperature  and  a 
Gaussian  distribution  in  the  logarithm  of  the 
relaxation  times  {lognormal  distribution)  proposed 
by  Yager  [22]  for  dielectric  relaxation  has  been  found 
by  most  investigators  to  fit  all  or  part  of  the 
ionic  diffusion  peak  in  several  glasses,  with 
activation  energies  around  30  kcals/mole  [16,  18,  19]. 
Nowick  and  Berry  [23]  have  worked  out  the  analysis 
of  continuous  relaxation  spectra  with  a Gaussian 
distribution  function  and  have  numerically  obtained 
and  tabulated  the  relationships  between  the  three 
parameters  of  a lognormal  distribution  function  and 
experimentally  measurable  functions.  Other 
relaxation  mechanisms  causing  a damping  peak, 
particularly  grain  boundary  viscosity,  in  metals 
and  crystalline  ceramics  have  been  studied  by 
Ke  [24,  25],  Wachtman  and  Lam  [26],  and  Chang  [27], 
among  others. 
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The  damping  in  porcelain  enamels  in  their 
softening  temperature  range  is  very  high,  and  has 
hitherto  not  been  studied  except  for  crude  experiments 
which  merely  confirmed  the  damping  potential  of 
enamel  coatings.  However,  various  other  properties 
of  porcelain  enamels,  like  strength,  resistance  to 
corrosion,  erosion  and  thermal  shock,  have  been 
studied  and  the  techniques  for  application  of  the 
coating  have  been  perfected  [28] . Cowan  et  al  [29] 
also  reported  the  variation  of  the  Young's  modulus 
of  porcelain  enamels  with  temperature  but  did  not 
make  any  detailed  observations  in  the  softening 
temperature  range. 

This  thesis  describes  a test  program  designed 
to  study  the  damping  properties  of  porcelain  enamel 
coatings  in  the  softening  temperature  range.  The 
design  of  three  pieces  of  apparatus  for  the 
measurement  of  the  complex  modulus  of  thin  coatings 
with  softening  temperatures  in  the  range  400  to  750  C 
is  discussed.  The  frequency  range  covered  is  1 to 
600  Hz.  The  methods  for  determining  both  the 
complex  Young's  modulus  E*  and  the  complex  shear 
modulus  G*  are  shown.  A Gaussian  distribution  in 
the  logarithm  of  relaxation  times  is  shown  to  be  an 
adequate  description  of  the  relaxation  process.  The 
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softening  process  is  shown  to  follow  an  Arrhenius 
relationship  with  anomalously  high  apparent 
activation  energies  which  can  only  be  explained  by 
extensive  cooperative  processes. 

Polymers  exhibit  a transition  range  of 
temperature  and  stress  relaxation  studies  in  this 
region  have  been  made  [30],  Bartenev  [31]  suggests 
that  inorganic  glasses  can  be  looked  upon  as  having 
a polymeric  structure;  that  they  can  be  viewed 
as  highly  cross-linked,  spatial  polymers.  However, 
certain  differences  in  the  structure  and  behavior 
of  polymers  and  that  of  porcelain  enamels  exist 
and  will  be  discussed  in  Chapter  2. 

This  thesis  uses  the  concept  of  a distribution 
of  relaxation  times  and  the  data  reduction  methods 
applicable  to  polymers  and  applies  them  to  a 
specific  inorganic  material,  namely,  porcelain 
enamels.  It  demonstrates  that  such  a treatment  is 
appropriate,  in  spite  of  the  fact  that  porcelain 
enamel  is  chemically  and  structurally  different  from 
polymeric  materials.  These  methods  lead  to  several 
conclusions,  hitherto  unknown,  about  porcelain  enamels 
as  a class  of  materials.  Chapter  2 describes  the 
porcelain  enamel,  its  composition,  method  of  coating, 
structure,  the  mathematical  model  used  to  represent 
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its  response  in  the  softening  temperature  range,  and 
the  calculation  of  experimentally  measurable 
functions  from  the  model — the  complex  modulus,  in 
this  case.  The  equipment  used  to  study  the  material 
properties  of  porcelain  enamels  used  as  coatings  on 
base  metals  is  described  in  Chapter  3.  The  basic 
observation  was  the  decay  rate  of  a freely  vibrating 
system  as  a function  of  frequency  and  temperature, 
from  which  the  complex  modulus  of  the  porcelain 
enamel  was  determined  as  a function  of  temperature 
and  frequency  using  the  equations  developed  in 
Chapter  3.  Chapter  4 presents  the  measured  data 
and  compares  the  experimentally  observed  values  of 
the  complex  moduli  with  those  predicted  from  the 
assumed  model  of  Chapter  2.  A summary  and  discussion 
of  the  results  is  given  in  Chapter  5. 

1.2  Damping  Units 

A number  of  different  parameters  have  been  used 
as  quantitative  measures  of  the  damping  capacity  of 
materials.  Some  of  the  more  important  of  these 
will  be  defined  in  this  section. 

Damping  is  usually  defined  in  terms  of  the  energy 
dissipated  per  cycle  when  a material  is  subjected  to 
cyclic  strain  varying  sinusoidally  in  time.  The 
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specific  damping  capacity  or  loss  coefficient  or  loss 
factor  r)  is  given  by 


n = (AU)/(2ttu)  (1.1) 

where  Au  = energy  dissipated  per  cycle 
and  U = maximum  strain  energy  stored. 

If  the  sinusoidal  strain  results  solely  from  a uni- 
axial stress, 

AU  = Pads  = £aedt 
and  U = (1/2) Ee l 
where  e = de/dt 

e0  = strain  amplitude 
E = Young's  modulus 

and  if  indicates  integral  over  a complete  cycle. 

If  AU  and  U are  computed  for  a system,  then  equation 

(1.1)  gives  the  system  loss  coefficient. 

The  energy  dissipation  causes  a decay  or  a 

reduction  in  the  amplitude  of  vibration  of  the  system 

which  leads  to  another  measure  of  damping,  viz,  the 

decay  rate.  This  is  expressed  in  terms  of  the 

logarithmic  decrement,  6,  given  by 

6 = lnfx0/x  ) /n  (1.2) 

u n 

where  xc  = amplitude  at  a certain  cycle 
and  x^  = amplitude  n cycles  later. 


A third  measure  of  damping  is  the  quality  factor 
Q,  determined  from  the  resonant  curve  of  amplitude 
versus  forcing  frequency;  it  is  given  by 

Q'1  = (Af)/fr  (1.3) 

where  Af  = width  of  the  response  curve  at 
an  amplitude  of  1//2  of  maximum 
and  fr  = resonant  frequency 

For  systems  with  small,  linear  damping, 
independent  of  amplitude,  it  can  be  shown  that  [1] 

n = Q'1  = 6/tt  (1.4) 

Oberst  [29]  characterized  the  damping  behavior 
of  viscoelastic  materials  subjected  to  sinusoidal 
strain-time  deformation  by  means  of  a complex 
elastic  modulus  E*,  where  E*  = E'  + iE"  = E'  (1+in) . 
E',  the  storage  modulus,  is  the  ratio  of  that  part 
of  the  stress  which  is  in  phase  with  the  strain,  to 
the  total  strain;  E",  the  loss  modulus,  is  the  ratio 
of  that  part  of  the  stress  that  leads  the  strain  by 
90°,  to  the  total  strain;  and  n is  the  loss  factor 
or  loss  coefficient  = E"/E' . E'  is  related  to  the 
maximum  stored  energy  while  E"  is  related  to  the 
energy  dissipated  per  cycle.  For  the  uniaxial  stress 
condition  above,  it  can  be  shown  that 


u = (1/2) E ' e* 

and  AU  = 7rE"e* 

0 

therefore,  (AU)/(2ttU)  = E"/E'  = n 
so  that  n defined  in  terms  of  the  real  and  imaginary 
parts  of  the  complex  modulus  is  the  same  as  the  n 
defined  for  a material  in  terms  of  the  energy 
dissipation  and  maximum  stored  energy  per  cycle. 


Chapter  2.  Porcelain  Enamels 


J 


2.1  Introduction  j 

Porcelain  enamels  are  thin  coatings  used  on  metal 
surfaces  for  decorative  and  protective  purposes. 

i 

Andrews  [6]  gives  a detailed  description  of  the 
enamels,  their  properties,  uses  and  application 
techniques.  A brief  physical  description  of  the 
enamels,  their  composition  and  method  of  coating  is 
given  in  Appendix  (1)  . 

2.2  Structure  of  Porcelain  Enamels 

There  is  little  information  available  in  the 
literature  regarding  the  microstructure  of  porcelain 
enamels.  They  possess  an  amorphous  structure  with 
little  order  in  it.  The  major  constituent,  silica, 
of  porcelain  enamels,  usually  present  up  to  about 
30%,  possesses  a vitreous  or  glassy  structure  with  only 
short  range  order.  Even  this  is  broken  up 
considerably  by  the  numerous  other  substances  that 

make  up  the  porcelain  enamel  (Appendix  (1)).  In  its  ; 

simplest  form,  porcelain  enamel  is  a highly  complex 
glass . 

Qualitatively,  it  is  a multi-component,  multi- 
phase system;  some  of  the  various  oxides  present  are 
in  solution  in  the  glassy  matrix  while  others  are 
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held  in  suspension.  The  suspended  phases  themselves 
are  often  crystalline.  The  components  of  the  enamel 
are  very  inert  and  do  not  react  with  one  another 
chemically.  The  matrix  usually  consists  of 
vitreous  silica.  The  oxides  held  in  suspension 
include  those  of  titanium,  tin,  antimony,  zinc, 
sodium,  potassium,  etc.  Other  chemicals  are  used 
in  the  enamel  to  ensure  that  the  suspended  particles 
are  evenly  dispersed  in  the  matrix.  An  important 
property  imparted  to  the  enamel  by  suspended  oxides 
like  those  of  tin  and  titanium,  is  opacity,  and  has 
been  the  object  of  chosen  scrutiny  [32], 

Porcelain  enamels,  because  of  their  vitreous 
nature,  do  not  possess  a definite  melting  point,  but 
are  characterized  by  a decreasing  viscosity  as  the 
temperature  is  increased.  The  temperature  band 
width  over  which  this  softening  takes  place  is 
relatively  wide.  In  this  temperature  range,  there 
are  large  changes  in  moduli  for  very  small  changes 
in  temperature,  and  it  is  in  this  temperature  range 
that  porcelain  enamels  are  useful  as  damping 
materials.  This  kind  of  transition  phenomena, 
characterized  by  rapid  modulus  changes,  has  been 
observed  and  studied  in  several  organic  polymeric 
materials  [33,34].  However,  porcelain  enamels, 
although  exhibiting  a transition  like  polymers,  are 
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very  different  materials,  both  chemically  and  physically. 
They  are  inorganic  as  opposed  to  the  organic  nature 
of  polymers.  Bartenev  [31]  has  suggested  that 
inorganic  glasses  do  indeed  possess  a polymeric 
structure — that  of  a highly  cross-linked  spatial 
polymer.  In  its  simplest  form,  porcelain  enamel  is 
a very  complex  glass,  with  only  as  much  as  30%  of  the 
glass  former,  silica.  Enamels  cannot  be  viewed  as 
multi-component  glasses.  It  is  true  that  they  are 
multi-component  systems,  but  they  are  also  multi- 
phase systems  which  is  not  the  case  for  multi- 
component  glasses.  Several  of  the  components 
present  in  porcelain  enamels  are  not  dissolved  in 
the  glassy  matrix;  they  are  merely  held  in 
suspension,  but  dispersed  uniformly.  For  example, 
the  opacity  of  porcelain  enamels  is  caused  by  the 
scattering  of  light  incident  on  crystals  of 
titanium,  tin  or  antimony  oxides  present  in  the 
enamels.  Porcelain  enamels  are  homogeneous 
mixtures  of  a large  number  of  inorganic  oxides  with 
varying  structures,  crystalline  or  amorphous;  the 
major  constituent  is  silica  which  has  a vitreous 
structure . 

Another  point  of  dissimilarity  between  porcelain 
enamels  and  organic  cross-linked  polymers  is  the 
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behavior  at  temperatures  greater  than  the  transition 


temperature.  Cross-linked  polymers,  in  general, 
show  a rubbery  region,  the  transition  being  from  a 
'glassy'  state  to  a 'rubbery'  state  [35] . Porcelain 
enamels,  on  the  other  hand,  transit  from  a 'glassy' 
state  to  a liquid  state;  they  behave  as  highly 
viscous  Newtonian  liquids  at  temperatures  greater 
than  the  transition  temperature.  This  behavior  is 
like  that  of  un-cross  linked  polymers,  which  also 
exhibit  a transition  from  a 'glassy*  to  a liquid 
state.  It  is  quite  apparent  from  their  other 
properties,  however,  that  structurally,  porcelain 
enamels  cannot  be  viewed  as  un-cross  linked  polymers. 

The  effect  of  the  numerous  oxides  suspended  in 
the  vitreous  silica  is  to  decrease  the  softening 
temperature  band  considerably.  Pure  vitreous  silica 
softens  in  the  range  1200  to  1500  C.  The  porcelain 
enamels  used  in  the  study  reported  in  this  thesis 
had  softening  temperatures  between  400  and  800  C. 
Table  2.1  shows  the  composition  of  two  enamel  frits 
used  in  the  tests.  The  difference  in  composition 
between  the  two  frits  is  small  - a 6%  lower  silica 
content  in  frit  #16845  is  made  up  by  slight  increases 
in  the  percentages  of  zinc  oxide,  magnesium  oxide, 
antimony  oxide  and  molybdenum  oxide.  However,  the 
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softening  temperature  bands  for  these  two  coatings 
differed  widely.  They  were  centered  around  465  C 
for  frit  #16845  and  around  600  C for  frit  #12325. 

As  one  would  expect,  the  softening  temperature  is 
higher  for  the  frit  with  the  greater  silica  content. 

Most  porcelain  enamels  have  specific  gravities 
varying  between  2.35  and  3.00.  Table  (2.1)  shows 
that  frit  #16845  contains  more  of  the  heavier 
elements,  zinc,  antimony  and  molybdenum,  than  does 
frit  #12325,  and  the  specific  gravities,  as  reported 
by  the  manufacturer  [36],  are  2.94  for  frit  #16845 
and  2.84  for  frit  #12325.  In  this  thesis,  it  will 
be  shown  that  most  porcelain  enamels  have  the  same 
value  for  elastic  moduli  below  their  softening 
temperatures . 

2.3  Structure  of  Glass 

As  was  mentioned  in  the  last  section,  the  major 
component  of  porcelain  enamels  has  a glass-like 
structure  and  it  is  this  component  which  is  responsible 
for  the  high  damping  capacity  of  porcelain  enamels  at 
elevated  temperatures.  Hence,  a brief  review  of  the 
structure  of  glass  in  general,  and  the  effect  of 
temperature  on  a glass,  is  in  order. 
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Glass  basically  consists  of  tetrahedral  units  of 
silica,  but  is  far  from  being  crystalline  [12].  Each 
unit  consists  of  an  atom  of  silicon  surrounded  by 
four  oxygen  atoms  to  form  a tetrahedron  with  the 
silicon  at  the  centroid  and  the  oxygen  atoms  at  the 
vertices.  Each  oxvgen  atom  is  shared  by  two 
neighboring  tetrahedral  units.  In  crystalline 
silicon  dioxide  (quartz) , the  relative  orientations 
of  these  basic  units  are  fixed.  In  glass,  however, 
this  orientation  is  random;  glass  possesses  short 
range  order  but  long  range  disorder.  The  structure 
of  glass  can  be  viewed  as  that  of  a liquid  "frozen  in. 

Some  authors  have  suggested  that  glass  should  be 
looked  upon  as  a fourth  state  of  matter  [37] , not  as 
a transparent,  fragile  substance,  but  as  a separate 
state  in  which  several  substances  can  be  made  to 
exist.  Some  substances,  like  silica,  lend  themselves 
readily  to  the  formation  of  a glass  and  are  called 
"glass  formers."  Others,  like  metals,  if  cooled 
rapidly  enough  from  the  liquid  state  and  if 
crystallization  can  be  prevented,  can  be  made  to  form 
a glass,  but  are  not  very  stable  in  this  form. 

2.3.1  Glass  Transition  Temperature 

To  understand  the  effect  of  temperature  on  a 
glass,  consider  how  a typical  property,  like  volume. 
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varies  with  temperature.  Figure  (2.1)  shows  this 
variation.  Starting  from  the  high  temperature  end 
of  the  figure,  on  cooling  a glass  forming  liquid, 
the  volume  at  first  decreases  like  any  other  liquid. 

At  the  melting  (freezing)  point,  shown  by  , most 
pure  crystalline  substances  show  a sudden  decrease 
in  volume.  The  volume  of  a glass,  however,  continues 
to  decrease  at  the  same  rate.  On  further  cooling, 
the  slope  decreases  rapidly  and  the  volume  continues 
to  decrease  at  a slower  rate,  which  is  often  the 
rate  at  which  the  solid  crystal  volume  would 
decrease.  This  behavior  is  sometimes  referred  to  as 
a "second  order  transition",  because  the  discontinuity 
is  not  in  the  volume  itself  but  in  its  derivative. 

The  temperature  range  over  which  this  transition 
takes  place  is  roughly  identified  by  the  "Glass 
Transition  temperature"  T . It  is  not  a well  defined 
temperature  like  the  melting  point  of  a crystalline 
solid  and  has  been  found  to  depend  on  the  rate  of 
cooling  [38]  . The  transition  temperature  is  lowest 
for  very  slow  cooling  rates  and  increases  with 
increasing  cooling  rates.  In  practice,  the  large 
times  involved  prohibit  the  measurement  of  very  low 
transition  temperatures.  It  is  conjectured  that 
there  is  a lower  limit  for  the  transition  temperature 
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and  that  it  is  10  to  100  C below  the  average 
experimental  value,  depending  upon  the  material. 
However,  the  transition  temperature  appears  to  increase 
without  any  limit  as  the  characteristic  time  is 
decreased.  Glass  transition  temperatures  well  above 
the  average  measured  have  been  measured  in  dynamic 
experiments  using  ultrasonic  frequencies  where  the 
characteristic  time  is  very  small  [39] . 

2.4  Relaxation  in  Porcelain  Enamels 

The  storage  Young's  modulus  E'  and  the  loss 
modulus  E"  for  a typical  porcelain  enamel  as  observed 
by  tests  described  in  Chapter  (3.1),  are  shown  as 
functions  of  temperature  in  figure  (2.2).  In  the 
transition  temperature  range  the  storage  modulus 
decreases  rapidly  with  increasing  temperature  and  the 
loss  modulus  attains  a maximum.  The  damping  in  this 
temperature  region  is  very  high. 

The  rapid  decrease  in  the  storage  modulus  is 
caused  by  the  relaxation  of  the  atomic  bonds  in  the 
vitreous  network  of  the  material.  Relaxation  is  the 
adjustment  of  a material  system  with  time  to  a new 
equilibrium  state  in  response  to  a change  in  an 
external  variable.  A characteristic  time  for  this 
adjustment  is  termed  the  relaxation  time.  A 
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definition  commonly  used  for  the  relaxation  time  is 
the  time  taken  for  the  modulus  of  the  material  to 
relax  through  1/e  of  the  total  relaxation.  The 
relaxation  time  decreases  with  increasing  temperature 
and  a transition  occurs  in  the  material  behavior 
when  its  temperature  has  been  sufficiently  increased 
so  that  the  relaxation  time  of  the  material  "crosses" 
the  time  scale  of  the  experiment  - the  period  of 
vibration  in  this  case. 

The  behavior  of  an  anelastic  material  is  often 
described  by  mechanical  models,  which  are  discussed 
in  Appendix  (2) . We  start  with  the  dependence  on 
frequency,  oo,  of  the  real  part  M'  and  the  imaginary 
part  M"  of  the  complex  modulus  M*  for  the  case  of  a 
system  with  a single  relaxation  time,  t.  The  real 
part  M'  decreases  from  My  to  MR  as  the  frequency 
decreases,  where  My  is  the  unrelaxed  modulus  and 
MR  is  the  relaxed  modulus.  The  imaginary  part  M" 
shows  a peak  at  a value  of  frequency  such  that  cot  = 1. 
When  M"  is  plotted  against  log  (cot)  , the  curve  has  a 
width  of  1.144  at  half  maximum.  This  is  called  a 

Debye  peak. 

The  experimentally  observed  peaks  in  the 
functions  E"  and  G"  for  porcelain  enamels,  when 
plotted  against  log  co,  are  much  broader  than  the  peak 
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appropriate  for  a system  with  a single  relaxation  time. 
This  suggests  a spectrum  of  more  than  one  relaxation 
time,  rather  than  a single  relaxation  time.  This  is 
also  reasonable  from  the  point  of  view  of  the  disorder 
in  the  structure  of  the  enamel.  A discrete  spectrum 
of  relaxation  times  can  be  obtained  by  choosing  a 
multi-parameter  (more  than  three)  mechanical  model  to 
represent  the  material.  But  this  is  not  a very  useful 
method,  as  there  is  no  indication  from  the  experimentally 
observed  modulus  functions  as  to  the  exact  values 
of  the  relaxation  times  or  their  corresponding 
relaxation  magnitudes.  The  functions  E"  and  G" 
could  be  used  to  determine  a discrete  spectrum  using 
Tobolsky's  procedure  X [30].  This  method  can  be  used 
to  obtain  meaningful  relaxation  times  only  if  the 
relaxation  times  are  widely  spaced,  at  least  a 
decade  apart.  Since  there  is  no  evidence  of  such 
being  the  case,  and  there  is  reason  to  believe  from 
the  structure  of  the  enamel,  that  a continuous 
relaxation  spectrum  would  be  appropriate,  we  shall 
assume  a continuous  spectrum  of  relaxation  times  for 
the  porcelain  enamel. 

For  an  anelastic  material  with  a distribution 
function  <J>(  &n  t)^  for  the  logarithm  of  the 

A logarithmic  scale  for  time  and  frequency  turns 
out  to  be  particularly  convenient. 


relaxation  times  x,  the  storage  modulus  M'  and  the  loss 
modulus  M",  from  [40],  are 

M'  (co)  = Md+6M/  (<{)(X.n  t)  u>2x  2 ) / (1+(jj2t  2 ) d£nx  (2.1) 

— 00 

00 

and  M"  (to)  = 6M/  (<M^n  x)  wx ) / (l+to2x 2 ) d£nx  (2.2) 

— 00 

where  SM  = the  magnitude  of  the  relaxation 
= ^ ~ MR 

and  $(in  x)  is  a distribution  function 

normalized  so  that 

CO 

I <p(ln  x)d£nx  = 1 (2.3) 

— 00 

Nowick  and  Berry  [40]  have  shown  that  the 
functions  M'  and  M"  are  not  particularly  sensitive 
to  the  form  of  the  distribution  function  $(£n  x). 

For  the  porcelain  enamels,  we  shall  assume  a lognormal 
distribution  of  relaxation  times,  i.e.  a Gaussian  or 
normal  distribution  in  the  logarithm  of  the  relaxation 
times.  Such  a distribution  is  completely  specified 
by  three  parameters,  namely,  the  magnitude  5m  of  the 
relaxation,  the  mean  relaxation  time  x , and  the 
width  B of  the  distribution.  It  is  a realistic 
distribution  and  at  the  same  time  mathematically 
tractable.  Nowick  and  Berry  have  worked  out  the 
mathematics  of  this  distribution  which  is  briefly 
shown  in  Appendix  (2) . Since  the  relaxed  modulus  MR 
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For  many  thermally  activated  relaxation  processes 
the  relaxation  times  have  been  found  to  follow  the 
Arrhenius  relation 

T = exp  (Q/RT)  (2.4) 

where  is  an  inverse  frequency  factor, 
independent  of  temperature. 

Q is  the  activation  energy  of  the  process, 
independent  of  temperature. 

R is  the  gas  constant 

= 1.986  cals./mole/°K 
and  T is  the  absolute  temperature. 

A distribution  in  the  relaxation  times  x can  be 
caused  either  by  a distribution  in  or  by  a 
distribution  in  Q,  or  both.  The  assumption  that  the 
distribution  is  only  in  the  parameter  is 
equivalent  to  assuming  that  the  distribution  peak 
width  parameter  8 is  independent  of  temperature.  This 
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is  a valid  assumption  because  the  temperature 
differences  that  are  required  to  change  the 
relaxation  time  by  several  decades  are  very  small. 
Typically,  for  porcelain  enamels,  a frequency 
increase  by  a factor  of  10  shifts  the  M"  function 
by  about  30F  (17C)  on  the  temperature  axis.  If  we 
assume  that  Q is  not  a distributed  function,  the 
mean  relaxation  time  is: 


Tm  = T*m  exP(Q/RT) 


(2.5) 


where  x.  is  the  mean  of  x. 


If  x . is  the  mean  relaxation  time  at  an 
mO 

arbitrarily  chosen  reference  temperature  Tn,  then: 


Tmn  = T*m  exP  (Q/RTn) 


(2.6) 


and , 


!°g  (*mAm0)  = (Q/2.303R)  (I/T-1/Tq) 


log  aT  = (Q/2.303R) (1/T-1/T0) 


(2.7) 


where  aT  is  the  ratio  of  the  mean  relaxation  time 
xm  at  temperature  T to  the  mean  relaxation  time  x g 
at  temperature  TQ.  Thus,  the  curve  of  log  aT  vs.  T 1 
is  a straight  line  whose  slope  is  (Q/2.303R). 

For  an  anelastic  material,  Zener  [41]  has  shown 


that  the  maximum  damping  occurs  when 


w • xm  = 1 
m m 

Then  uim  is  the  frequency  in  rad/sec  for  maximum 
damping  at  temperature  T,  so  that 
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(2.8) 


O'  = T /T  a 

T m mO 


= V%  = fO/fro 


where  fQ  is  the  frequency  in  Hz  for  maximum  damping 
at  the  reference  temperature  Tq,  and  is  the 
frequency  in  Hz  for  maximum  damping  at  temperature 
T.  Equations  (2.7)  and  (2.8)  provide  a means  for 
determining  the  activation  energy  Q of  the  relaxation 
process  from  experimentally  observed  M’  or  M"  data 
as  functions  of  temperature  at  various  frequencies. 
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Chapter  3.  Experimental  Technique  and  Calculations 


3.1  Experimental  Technique 

A number  of  different  parameters  have  been  used 
as  quantitative  measures  of  damping,  as  indicated  in 
Chapter  (1.2).  The  one  most  often  used  for 
engineering  purposes  is  the  energy  dissipated  per 
cycle  when  a material  is  subjected  to  sinusoidal 
strain-time  deformation.  The  energy  dissipation 
leads  to  a reduction  in  the  amplitude  or  decay  of 
vibration  of  the  system.  The  most  commonly  used 
method  for  determining  the  damping  capacity  of  a 
system  is  the  measurement  of  the  rate  at  which  the 
amplitude  of  free  vibration  of  the  system  decreases . 
Damping  capacity  can  also  be  determined  from  the 
increase  in  amplitude  of  a vibrating  system  as  the 
excitation  frequency  is  varied  through  resonance 
holding  the  amplitude  of  the  exciting  force  constant. 
The  energy  dissipation  can  also  be  determined  by 
measuring  the  power  necessary  to  maintain  vibration 
at  a given  amplitude. 

The  decay  rate  of  a freely  vibrating  system  was 
the  basic  observation  in  the  method  described  in 
this  thesis.  Since  it  is  very  difficult  to  obtain  a 
wide  freqeuncy  range  from  any  one  simple  test 
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configuration,  three  different  pieces  of  apparatus 
were  designed  to  cover  the  frequency  range  from 
1 Hz  to  600  Hz.  They  were  designed  to  operate 
satisfactorily,  without  undergoing  permanent  damage, 
in  the  temperature  range  70F  (20C)  to  2000F  (1100C) . 
A torsion  pendulum  specimen  was  used  for  frequencies 
between  1 Hz  and  15  Hz.  A cantilever  beam  specimen 
was  used  in  the  range  10  Hz  to  170  Hz,  and  a hollow 
cylindrical  torsion  specimen  was  used  in  the 
frequency  range  165  Hz  to  600  Hz.  In  each  case,  the 
frequency  of  free  vibration  was  changed  by  changing 
the  total  mass  of  the  system.  The  specimens  were 
made  of  high  temperature  alloys  and  were  coated  with 
the  porcelain  enamel  over  part  of  the  surface. 
Measurements  were  made  on  each  configuration  at  a 
number  of  different  temperatures  in  the  range  of 
interest.  All  parts  of  the  apparatus  subjected  to 
high  temperatures  were  made  of  the  same  material  in 
order  to  avoid  geometry  changes  due  to  differences 
in  thermal  expansion.  Heating  was  by  electric 
furnaces  and  the  coated  sections  of  the  specimens 
were  maintained  at  a uniform  temperature.  In  the 
v ^st  case,  the  temperature  variation  across  the 
coating  was  less  than  20F  (11C) . Since  system 
measurements  depend  upon  integrated  properties  of 
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the  whole  volume  being  tested,  the  apparent  dependence 
on  temperature  is  reduced  if  the  various  parts  of  the 
specimen  are  at  different  temperatures.  The  effective 
bandwidth  in  the  present  tests  was  about  100F  (55C) , 
and  20F  (11C)  is  about  the  maximum  variation  that 
can  be  tolerated. 

In  each  case,  a transducer  was  used  to  sense  the 
mechanical  vibration  and  convert  it  to  an  analogous 
voltage.  This  voltage  was  used  to  drive  a recorder 
to  determine  the  decay  rate,  and  a frequency  counter 
for  the  frequency.  The  details  of  the  design  of  the 
apparatuses,  the  instrumentation  used  and  the 
calculation  of  the  natural  frequencies,  for  the 
three  cases  are  discussed  in  Appendices  (3),  (4)  and 

(5)  respectively.  Appendix  (6)  describes  the  various 
procedural  problems  that  were  encountered  during  the 
tests  and  the  measures  taken  to  overcome  them. 

3.2  Calculations 

The  complex  elastic  modulus  of  the  coating 
material  is  found  from  the  measured  system  damping, 
frequency  and  specimen  geometry  in  all  three  systems. 
The  general  method  is  the  same,  but  the  details  are 
somewhat  different  in  each  case.  Subscripts  b and  c 
will  be  used  to  denote  system  quantities;  b for  the 
system  with  a bare  specimen  and  c for  the  system  with 
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a coated  specimen.  Subscripts  1 and  2 will  be  used 
to  denote  the  materials;  1 for  the  base  material 
and  2 for  the  coating  material. 

The  complex  Young's  modulus  E*  of  the  coating 
material  is  determined  from  the  bending  test  measure- 
ments following  the  methods  used  by  Van  Oort  [42] 
and  Nashif  [43] . The  following  assumptions  are  made 
in  deriving  the  formulae  for  determining  E*. 

1)  The  beam  is  of  uniform  thickness  as  shown 
in  figure  (3.1) . 

2)  Anticlastic  curvature  is  completely 
constrained  so  that  s =0. 

ZZ 

3)  Since  the  coating  is  very  thin,  and  the 
normal  stress  on  the  surface  is  zero, 

cr  is  assumed  to  be  zero  in  the  coating. 

4)  The  base  metal  damping  is  zero,  so  that 
the  Young's  modulus  of  the  base  metal 
is  real. 

5)  The  kinetic  energy  of  the  mass  at  the  free 
end  is  far  greater  than  the  kinetic  energy 
of  the  beam  itself.  The  mass  of  the 
coating  itself  is  so  small  that  it  does 
not  appreciably  change  the  kinetic  energy 
of  the  system. 

6)  All  strain  energy  stored  is  in  the  beam 
itself  and  none  in  any  other  part  of  the 
system. 
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Following  the  same  development  as  Lee  [8] , the 


complex  bending  moment  M*  of  the  system  from  which 
the  system  loss  coefficient  can  be  determined  is 


M*  = M'  + iM" 
c c c 

= (2/3p) [EjdJ+E£ (3d^d2+3d1d2+d2> ] 
+ inEE^(2/3p)  (3djd2+3d1d2+d2) 


where 


p = radius  of  curvature  of  the  beam 


d^  = half  thickness  of  the  specimen 
62  = thickness  of  the  coating 
E^  = Young's  modulus  of  base  metal 
EJ,  = Real  part  of  Young's  modulus  E*  of  the 
coating  material 

nE  = loss  tangent  of  the  coating  material 

- E2/E2 

E2  = Imaginary  part  of  Young's  modulus  E* 
of  the  coating  material. 

The  loss  coefficient  for  uniform  bending  in 
this  system  is 

n = m;/m* 

= nE  ac/ (1+ac)  (3.1) 

where 

a = modulus  ratio  E2/E^ 
c = 3t  + 3t2  + t3 
t = thickness  ratio  d2/d^ 


The  frequency  of  the  system  with  a bare 
specimen  is  given  by 

fb  - V<KE) 

where  = the  strain  energy  stored  and  (KE)  is  a 
kinetic  energy  term.  For  the  cantilever  specimen 
with  a large  rigidly  connected  mass  at  the  free  end, 
the  kinetic  energy  of  the  mass  is  much  greater  than 
that  of  the  beam  itself. 

Now, 

U,  = / (1/2)0  e dv 
b v xx  xx 

- 2E"  /o1(1''2IElExx 
2 ^12 

= (LwE^/ p ) f y ay 


= (L/2p":)E1I1 


where 


L = length  of  the  beam 
w = width  of  the  beam 
p = radius  of  curvature 
and  1^  = area  moment  of  inertia  of  the  cross 
section  of  the  bare  beam 
= ( 2/3 ) d^w 
Therefore , 


fb  = “l1! 


where 


K = L/ ( 2 p ( KE ) ) 


For  the  system  with  the  coated  specimen. 


fc  - V(KE) 


Again, 


U = f (l/2)o  e dv 

c v ' xx  xx 

d +d 

= 2Lw  / (1/2)0  e dy 

Q XX  XX  1 

di  2 "*"^2  2 

= LW  ;Q  Elexxd*  + Lw  E2£xxd^ 

= ( L/ 2 p 2 ) (E1I1  + E'I2) 


where 


I2  = area  moment  of  inertia  of  the  coating 


= (2d^d2  + 2d1d2  + (2/3)d2)w. 

Therefore , 

£c  - K<E1I1  + E2J2> 

where  K is  the  same  constant  as  for  the  uncoated 
specimen.  Hence, 

EX  * 1 * E2V(EiV 

= 1 + al2/lx  (■ 

I2/I1  = (2djd2  + 2dxd^  + (2/3)d2/(2/3)dJ 


= 3d2/d1  + 3d*/d*  + dd/d^ 

2 3 

* 3t  + 3t  + tJ 


Equation  (3.2)  becomes 


2 2 

f /ff  - 1 + ac 
c'  b 

Substituting  into  equation  (3.1) 
nE  = n(l  - f l/ty1 


(3.3) 


(3.4) 


Equation  (3.2)  gives 


E2  ' <VC>  <£c/£b  - 11 


(3.5) 


and  from  the  definition  of  n. 


E2  ‘ E2nE 

* <VC)  n£c/fb 


(3.6) 


Equations  (3.5)  and  (3.6)  are  valid  for  uniaxial 

stress  conditions  but  our  bending  specimen  has  a 

length  to  width  ratio  of  almost  one;  the  anti- 

clastic  curvature  can  be  assumed  to  be  fully 

constrained  and  a Poisson's  ratio  correction  of 
2 

(1-v  ) is  necessary. 

i.e.  E'2  = (1-v2)  (Ex/c)  (f2/f2-l)  C 


E^  = (1-v2)  (Ej/c) nf^/f£ 


(3.8) 


The  complex  shear  modulus  G*  of  the  coating 
material  is  determined  from  the  torsion  pendulum 
tests.  The  following  assumptions  are  made  in  this  case, 
1)  The  coating  material  is  in  pure  shear  and 

the  shear  strain  amplitude  in  the  coating  is 
uniform. 
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2)  The  base  metal  damping  is  zero  so  that  the 


shear  modulus  of  the  base  metal  is  real. 

3)  The  kinetic  energy  of  the  specimen  itself 
is  negligible  when  compared  with  the 
kinetic  energy  of  the  mass  at  the  free 
end,  and  the  mass  of  the  coating  is  so 
small  that  it  does  not  change  the  kinetic 
energy  of  the  system. 

4)  All  strain  energy  is  in  the  length  of  the 
specimen  between  the  support  and  the  weight 
at  the  free  end. 

Equation  (3.5)  modified  by  a length  correction 

gives 

G2  = (Gj/c) (l-f£/fc>/(b-(l-fb/fc))  (3.9) 

where 

= shear  modulus  of  base  material 
c = 3t  + 3t^  + t3  where  t = 
r 2 = coating  thickness 

= radius  of  bare  specimen 
b = length  ratio 

= (coated  length  of  specimen) / (total  length). 
The  loss  modulus  is  determined  from  the  energy 
dissipated  per  cycle.  The  system  loss  coefficient  n 
is  defined  as 
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n = (AU)s/(2ttUs) 
where 

(AU)g  = energy  dissipated  per  cycle  and 
Ug  = the  maximum  energy  stored. 

For  the  torsion  specimens,  the  maximum  energy  stored 
is  equal  to  the  maximum  kinetic  energy,  which  is 
assumed  to  be  completely  in  the  mass  at  the  free 
end.  Thus 

Ug  = (1/2)  u)2j92  (3.10) 

where 

io  = circular  frequency 
J = mass  moment  of  inertia  of  the  mass 
9 = angular  displacement  of  the  mass 
The  strain  amplitude  in  the  coating  is  uniform 
through  its  thickness  and  is 

Y 0 - k0  (3.11) 

The  strain  was  measured  in  the  high  frequency  torsion 

specimen  with  strain  gages  and  the  constant  k was 

-2 

found  to  be  3.0  x 10  (figure  (A5.4)).  k was 
assumed  to  be  the  radius  to  length  ratio  3.3  * 10  ^ 
for  the  torsion  pendulum  specimen. 

If  there  is  no  damping  in  the  base  material, 
the  energy  dissipated  per  cycle  is, 

(AU)  g = ttG^VYq  (3.12) 

where 

'J  = volume  of  the  coating. 
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= (nJw2)/Vk2  (3.13) 


The  system  loss  coefficient  was  determined  from 
decay  measurements.  For  a linear  system,  with  low, 
amplitude  independent  damping,  in  free  vibration  [1] , 
ns  = (in(xQ/xn)  )/(irn) 
where 

Xq  = amplitude  at  a certain  cycle 
xn  = amplitude  n cycles  later. 

Now, 

*n(x0/xn)  = (in  10)  A(dB)/20 
= 2.303  A/20 

and 

n = ft 
where 

f = natural  frequency 
and 

t = time  for  a decrease  of  A dB. 

Therefore , 

ns  = ( 2 . 303/  (20tt  ) ) A/ft 
= 3.66  x 10~2  A/ft 

In  deriving  the  formulas  for  the  complex  moduli, 

it  was  assumed  that  energy  is  dissipated  entirely  in 

the  coating  and  that  there  is  no  damping  in  the  base 

material.  i.e.  ri  in  equations  (3.8)  and  (3.13)  is 

(AU)  / (2nU  ) where  (AU)  is  the  energy  dissipated  per 
OS  c 
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cycle  in  the  coating  and  Ug  is  the  maximum  energy 
stored  in  the  system.  This  is  not  strictly  true 
and  a correction  must  be  made  for  the  actual  system 
damping  with  a bare  specimen. 

Let 


nb  = damping  coefficient  for  the  system  with 
a bare  specimen 
= (&U)b/(27iUb) 

_ energy  dissipated/cvcle  in  bare  specimen 
2tt  (maximum  energy  stored) 

ns  = system  damping  coefficient  with  a coated 

specimen 

= (AU)  / (2nU  ) 

b b 

and 

n = (AU)  /(2ttU  ) is  the  damping  that  would  have 
c c s 

been  observed  if  the  bare  specimen  damping 
were  zero. 

Since  the  energy  stored  in  the  coating  is  negligibly 
small,  we  assume  that  Ub  = Ug,  and  that  (AU)c  = 

(AU) s - (AU)b.  Then, 

nc  = ((Au)g  - (AU)b)/(27iUs) 

= (AU)s/(2ttUs)  - (AU)b/(2TTUb) 

= ns  ' nb 

The  determination  of  the  complex  modulus  of  the 
coating  material  depends  critically  upon  the  ratio  of 
the  frequencies  of  the  coated  and  uncoated  specimens. 


36 


This  can  be  seen  in  equations  (3.7)  and  (3.9)  which 
use  the  small  differences  in  the  frequencies  of  the 
coated  and  uncoated  specimens.  In  some  cases,  the 
bare  specimen  was  tested  before  coating  it.  Since 

there  was  little  difference  among  the  measured  ! 

values  of  damping  coefficients,  an  average  curve  of 

loss  coefficient  as  a function  of  temperature  was 

assumed  for  all  specimens,  except  for  support  effects 

which  were  assumed  to  be  independent  of  temperature. 

The  frequency  change  with  temperature  was  assumed  to 
depend  only  on  the  modulus  change  and  thus  to  be  in 
the  same  proportion  for  all  specimens. 

At  a sufficiently  high  temperature  the  modulus 
of  the  coating  decreases  enough  that  its  contribution 
to  the  stiffness  of  the  specimen  is  negligible.  For 
thin  coatings,  then,  the  frequency  of  the  system 
with  the  coated  specimen  is  the  same  as  that  with 
the  bare  specimen  alone.  Hence,  the  frequency  curve 
of  the  bare  specimen  is  adjusted  proportionally 
so  that  the  frequency  with  the  coated  specimen 
approaches  that  with  the  bare  specimen  near  the 
firing  temperature  of  the  coating.  This  adjusted 
curve  is  used  as  an  estimate  of  the  bare  specimen 
frequencies  of  the  given  specimen. 
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It  is  assumed  that  there  is  no  energy  dissipation 
in  the  coating  at  room  temperature  and  that  any  change 
in  the  measured  system  loss  coefficient  at  room 
temperature  between  the  coated  and  uncoated  specimens 
is  due  to  support  effects  and  is  independent  of 
temperature.  The  bare  specimen  damping  curve  is 
adjusted  vertically  to  match  the  curve  for  the 
coated  specimen  at  room  temperature.  The  bare 
specimen  loss  coefficient  for  a given  specimen  is 
found  from  this  curve.  Figure  (3.2)  shows  the  system 
loss  coefficient  and  frequency  as  functions  of 
temperature  for  a typical  specimen.  The  adjusted 
bare  specimen  damping  and  frequency  curves  are  also 
shown . 

In  the  bending  tests  at  11  Hz  and  47  Hz,  the 
heavy  arm  rigidly  attached  to  the  free  end  of  the 
specimen  justifies  the  assumption  that  strain 
energy  is  stored  only  in  the  beam  itself  and  none 
in  the  heavy  mass.  However,  in  order  to  make 
measurements  at  170  Hz,  the  heavy  arm  was  replaced 
by  a light  rod  and  the  above  assumption  is  no 
longer  valid.  Figure  (3.3)  shows  the  loss 
coefficient  and  frequency  as  functions  of  temperature 
for  the  same  specimen  as  in  figure  (3.2),  when 
decaying  freely  at  170  Hz.  The  maximum  loss 
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coefficient  is  6.0  x 10  at  170  Hz  as  compared  with 
11.5  x 10  ^ at  47  Hz.  The  percentage  decrease  in 
frequency  as  a result  of  the  softening  of  the  coating 
is  less  for  the  high  frequency  measurement  than  it  is 
for  the  lower  frequency  one.  If  we  assume  that  strain 
energy  is  stored  in  the  specimen  and  elsewhere  in  the 
system  also  at  170  Hz,  but  only  in  the  specimen  at 
11  or  47  Hz,  then  we  may  correct  the  damping 
coefficient  to  get  comparable  values  for  the  170  Hz 


case . 


Let  U 1 be  the  energy  stored  in  the  specimen 


alone 


u(1)  = + u2(1) 

where 


U1(l)  = energy  stored  in  the  base  metal 


U2  ' ' - energy  stored  in  the  coating. 


Let  IP  ' be  the  energy  stored  elsewhere  in  the 

(2) 

system.  Then  for  a lower  frequency  fA,  UA  =0, 
where  the  subscript  A refers  to  a lower  frequency, 
11  Hz  or  47  Hz. 

Let  f be  the  frequency  that  would  have  been 
observed  if  the  coating  had  not  softened  rapidly. 


Then 


f A " '“li1’  + U2A1,,/itE(C0) 
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where  KE(Cq)  is  a function  of  the  strain  amplitude 
e„  only. 

Let  f be  the  frequency  actually  measured  when 
A 

the  coating  softens.  Then 

~fl  - uJ1,/KE(Eo> 

Therefore , 


2 ~ 2 ~ 2 

r a $ (fA  - 


A'  ' A 


Now, 


“ 2(fA  ' fA>/fA 

= n (1) /u  (1)  0.14) 

U2A  /U1A 


n = loss  coefficient  at  a lower  frequency 
A 

where 

E"  * the  loss  modulus 
2A 

and 

V = volume  of  coating 


Then 


Ek  - raiAl\ 

where 


(3.15) 


k = 2/(eJv) 

At  the  higher  frequency  f0  = 170  Hz, 

U O)  = energy  stored  in  the  system  other  than 
B 

in  the  specimen 

* 0 


I 
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Defining  fQ  and  f0  similarly  to  fA  and  fA,  we  have 
fB  - '"IB1’  + "2B1’  + UB2),/KE(e0I 


*B  ‘ '“ill1’  + "I2'  >/KE(c0) 


Therefore , 


,,2  ~2,  ,t2 

{fB  ‘ V/fB 


2^fB  _ fB,/,fB 


" U2i1)/(U1B1)  + Ui2)) 


(3.16) 


nD  = loss  coefficient  at  frequency  fD 

D ID 


= Kn  / k(U1^1)  + U<2) ) 


Therefore , 

E2B  - "'"ll1’  + 

Since  is  independent  of  frequency 


(3.17) 


0 (1)  - U (D  = u(1)  1 

U1A  U1B  1 


u (i)  = u <D  = u(1) 

U2A  U2B  U2 


for  the  same  e. 


(2)  (2) 

and  writing  Ug  as  U , we  have  from  equations  (3.1) 


rA/rB  = R 


= (u;[1)  + u (2)  )/u|1) 


= 1+u^Vuj1^ 


= (l  + u(2)/u{1})  • nB/nA 
= r • nB/nA  (3.19) 

The  ratio  n^/n,  can  be  determined  from  the  measured 
B A 

values  of  the  loss  coefficient. 

Table  (3.1)  lists  the  frequencies,  frequency 
changes,  loss  coefficients  and  the  ratios  R and 
nA/nB  for  several  specimens  coated  with  porcelain 
enamels  of  various  types.  Since  the  ratios  R and 
r)A/oB  are  both  approximately  1.8  for  all  specimens, 
equation  (3.19)  shows  that  the  loss  modulus  of 
the  coating  material  does  not  depend  appreciably 
upon  frequency.  Therefore,  the  lower  value  of  the 
loss  coefficient  measured  at  the  higher  frequency  is 
not  due  to  an  increase  in  frequency  but  to  the 
change  in  the  test  conditions.  Hence,  the  loss 
coefficient  measured  at  170  Hz  should  be  multiplied 
by  1.8  to  obtain  the  loss  coefficient  that  would 
have  been  observed  had  all  of  the  strain  energy  been 
in  the  specimen  alone  and  none  in  the  rest  of  the 
system. 


Chapter  4.  Test  Results 


4.1  Introduction 

The  data  from  the  damping  experiments  conducted 
on  porcelain  enamel  coatings  using  the  method  and 
equipment  described  in  Chapter  3,  and  its  analysis 
will  be  presented  in  this  chapter.  The  system  loss 
coefficient  in  the  torsion  pendulum  and  bending  tests 
was  independent  of  stress  level.  The  torsion 
pendulum  tests  were  conducted  at  a shear  stress 
amplitude  of  5,000  psi,  and  the  bending  tests  at  an 
amplitude  of  10,000  psi.  For  the  high  frequency 
torsion  specimen  as  indicated  in  Appendix  (5) , the 
system  loss  coefficient  was  amplitude  dependent. 

The  results  of  tests  on  this  apparatus  correspond  to 
a shear  stress  level  of  50  psi.  However,  the  loss 
coefficient  in  the  coating  alone  did  not  depend  on 
the  stress  level.  The  various  coatings  tested, 
their  firing  temperatures  and  manufacturers  are 
listed  in  Table  (4.1).  Only  one  coated  high  frequency 
torsion  specimen  was  tested,  using  a frit  of  C.V. 
#12325. 

4.2  Presentation  of  Data 

Figures  (4.1)  to  (4.3)  show  the  measured  loss 
coefficient  and  frequency  as  functions  of  temperature 
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for  a torsion  pendulum  specimen  (#2)  at  three  nominal 
frequencies  of  1.2  Hz,  8.4  Hz  and  14.25  Hz.  The 
same  functions  are  shown  in  figures  (4.4)  to  (4.6) 
for  a bending  specimen  (#16A)  at  frequencies  of 
11  Hz,  47  Hz  and  170  Hz.  The  bare  specimen  damping 
and  frequency  curves  are  also  shown  in  each  case. 

Figure  (4.7)  shows  the  damping  nc  due  to  the  coating 
alone  and  frequency  as  functions  of  temperature  for 
the  high  frequency  torsion  specimen  at  a frequency 
of  165  Hz.  The  bare  specimen  damping  measured  in 
this  case  was  very  high,  greater  than  0.002  at  room 
temperature.  This  was  a result  of  improper  fit  in 
the  expansion  joints  between  the  test  section  and  the 
end  grips  of  the  specimen.  Tables  (4.2),  (4.3)  and 
(4.4)  list  the  maximum  loss  coefficient  and  correspond- 
ing temperatures  for  specimens  tested  in  low  frequency 
torsion,  bending  and  high  frequency  torsion 
respectively.  The  temperature  of  maximum  loss 
coefficient  was  determined  from  the  graphs  of  nc  vs 
temperature  for  all  specimens.  The  peaks  in  these 
functions  were  symmetrical  and  the  temperature  of 
maximum  loss  coefficient  was  obtained  as  the  average 
of  the  two  temperature  readings  corresponding  to 
the  same  value  of  nc  a little  below  the  maximum.  The 
tests  at  33  Hz  were  performed  by  Lee  [8]  and  the 
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results  shown  at  that  frequency  have  been  taken  from 
his  thesis.  The  temperature  points  in  Lee's 
thesis  around  the  damping  peak  were  very  widely 
spaced  and  the  maximum  loss  coefficient  and 
temperature  values  shown  are  not  very  reliable.  The 
coating  for  specimens  #16A  and  22,  bending,  used  a 
frit  with  the  same  number  (#12325)  as  that  for 
specimens  #16  and  18A,  bending;  but  was  prepared  and 
applied  at  a different  place  under  different 
conditions.  This  perhaps  explains  the  difference  in 
the  results  for  specimens  #16A  and  22,  when  compared 
with  those  for  specimens  #16  and  18A.  Specimens 
#15  and  15A,  bending,  were  coated  with  a frit 
carrying  the  same  number  (#16845)  also,  prepared  and 
applied  at  different  places;  but  the  results  in 
this  case  are  not  much  different.  Figures  (4.8)  and 
(4.9)  show  the  system  loss  coefficient,  nc  due  to 
the  coating  alone  as  a function  of  temperature  for 
two  torsion  pendulum  specimens  in  the  region  of 
maximum  damping  at  three  nominal  frequencies  1.2  Hz, 
8.4  Hz  and  14.25  Hz.  Figure  (4.10)  shows  the  same 
curves  for  a bending  specimen  (#16A)  at  freqeuncies 
11  Hz,  47  Hz  and  170  Hz.  The  loss  coefficient 
measured  at  170  Hz  has  been  multiplied  by  1.8,  as 
indicated  earlier,  to  obtain  the  loss  coefficient 
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that  would  have  been  measured  had  all  of  the  strain 
energy  been  in  the  bending  specimen.  Figure  (4.11) 
shows  nc  vs  temperature  for  the  high  frequency 
torsion  specimen  at  frequencies  165  Hz,  340  Hz, 

430  Hz  and  600  Hz.  Where  measurements  were  made  at 
sufficiently  high  temperatures  that  the  damping 
peak  could  be  completed,  the  peak  width  at  half 
maximum  loss  coefficient  is  also  shown  in  Tables 
(4.2)  to  (4.4) . 

The  following  general  observation  can  be  made: 

a)  All  coatings  exhibit  a transition  range  of 
temperature  during  which  the  system  loss  coefficient 
reaches  a maximum  and  the  frequency  drops  sharply. 
The  peak  in  the  loss  coefficient  occurs  at  a 
temperature  100  to  200C  below  the  firing  temperature 
of  the  coating. 

b)  For  a given  coating,  the  maximum  magnitude  of 
the  system  loss  coefficient  is  almost  directly 
proportional  to  the  coating  thickness,  at  least  for 
the  thicknesses  tested.  Figure  (4.12)  shows  n 

c max 

per  unit  thickness  vs  thickness  of  the  coating  for 
those  tested  in  bending. 

c)  For  coatings  made  from  different  frits  with 
approximately  the  same  thickness,  there  is  some 
indication  that  the  maximum  system  loss  coefficient 
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increases  with  the  maximum  damping  temperature.  For 
the  coatings  tested  in  bending,  figure  (4.13)  shows 
n per  unit  thickness  vs.  maximum  damping  temper- 

ature . 

d)  The  magnitude  of  the  maximum  system  loss 
coefficient  is  independent  of  frequency. 

e)  In  almost  every  case,  there  is  an  increase  in 
the  temperature  of  maximum  damping  with  increasing 
frequency.  This  shift  in  temperature  is  very  small, 
typically  about  20C  for  a frequency  factor  of  10; 
and  in  some  cases  this  temperature  shift  is  not 
apparent,  probably  because  of  small  errors  in  the 
temperature  measurement. 

f)  In  most  cases,  the  peak  width  at  half  maximum 
is  about  60C  in  bending  and  about  50C  in  the  torsion 
pendulum  tests. 


4.3  Data  Reduction  and  Analysis 


4.3.1  Complex  Modulus 


The  complex  modulus  of  the  coating  was  determined 
from  the  equations  derived  in  Chapter  (3.2).  The 
storage  shear  modulus  and  the  loss  shear  modulus 
of  the  coating  were  calculated  from  the  torsion 
pendulum  tests  using  equations  (3.9)  and  (3.13).  The 
storage  Young's  modulus  and  the  loss  Young's 
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modulus  were  calculated  from  the  bending  tests 
using  equations  (3.5)  and  (3.6).  The  results  from 
the  hollow  torsion  cylinder  apparatus  were  not 
satisfactory.  The  loss  shear  modulus  was 
calculated  for  coating  #12325  from  the  data  shown 
in  figure  (4.11)  obtained  using  the  hollow  torsion 
cylinder  apparatus.  The  value  obtained  for  G 
however,  was  inconsistent  with  the  loss  modulus 
values  from  the  bending  and  torsion  pendulum  tests 
for  the  same  coating.  Figures  (4.14)  and  (4.15) 
show  G^  and  °f  coating  #12325  as  functions  of 
temperature  calculated  from  tests  on  torsion 
pendulum  specimen  #2.  Figures  (4.16)  and  (4.17) 

I 11 

show  E 2 and  E^  for  the  same  coating  as  functions  of 
temperature  calculated  from  tests  on  bending 
specimen  #16A.  The  storage  modulus  shows  a sharp 
drop  with  increasing  temperature  in  the  transition 
region  and  the  loss  modulus  reaches  a maximum. 
Tables  (4.5)  and  (4.6)  list  the  room  temperature 
values  of  the  storage  modulus  and  the  maximum  loss 
modulus  values  for  the  various  coatings.  For 
torsion  pendulum  specimen  #3,  the  coating  was  so 
thin  that  the  system  frequency  difference  between 
the  coated  and  uncoated  specimens  was  very  small. 
Hence,  reliable  G^  values  for  coating  #16845  could 
not  be  obtained.  Except  for  results  from  specimen 
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#16,  the  Young's  storage  modulus  for  all  coatings  is 
about  10-11  =<  10^  psi  at  room  temperature.  The 
maximum  value  of  is  about  1.45  x 10^  psi  for 
coatings  #12325  and  #16845,  and  as  high  as  2.0  x 10^ 
psi  for  coating  S9919D.  For  coating  #12325,  the 
storage  modulus  and  the  loss  modulus,  obtained  from 
tests  on  sp.  #16,  are  both  about  30%  higher  than 
the  moduli  obtained  from  tests  on  sps.  #16A  or  22. 

The  coatings  on  the  three  specimens  carried  the 
same  frit  number,  but  the  frit  used  for  sp.  #16 
belonged  to  a different  batch  and  the  coating  on 
sp.  #16  was  applied  at  a different  place.  However, 
it  is  doubtful  that  this  is  the  cause  of  the  high 
modulus.  The  storage  modulus  of  other  porcelain 
enamels  tested  agree  with  the  values  obtained  from 
sps.  #16A  and  22.  The  high  values  obtained  for 
sp.  #16  are  probably  due  to  an  error  in  measurement. 
Equations  (3.5)  and  (3.6)  show  that  only  an  error 
in  the  measurement  of  the  coating  thickness  could 
cause  both  the  storage  modulus  and  the  loss 
modulus  to  increase  or  decrease  in  the  same 
proportion.  The  thickness  of  the  coating  on  sp.  #16 
is  only  0.00185  inch,  and  more  likely  to  be  in 
error  than  the  thicknesses  of  other  thicker  coatings. 
The  coating  thicknesses  were  obtained  as  the 
difference  in  thickness  between  the  coated  and 


uncoated  specimens.  Unfortunately,  once  a specimen 
is  coated,  the  uncoated  specimen  thickness  measure- 
ment cannot  be  repeated  to  check  the  value  obtained 
for  the  coating  thickness.  It  is  believed  that 
the  high  modulus  values  measured  for  the  coating 
on  sp.  #16  is  due  to  a 30%  error  in  the  measurement 

t 

of  the  coating  thickness.  Some  general  observations 
can  be  made  about  the  modulus  values. 

a)  The  maximum  value  of  the  loss  modulus  does  not 
depend  on  frequency. 

b)  The  width  of  the  loss  modulus  peak  at  half 
maximum  is  about  60C  for  coatings  tested  in  bending 

► 

and  50C  for  those  tested  in  torsion — the  same  as 
observed  for  the  nc  peak  in  the  two  cases. 

c)  There  is  an  increase  in  the  temperature  of 
maximum  loss  modulus  with  increasing  frequency.  This 
shift  in  temperature  with  frequency  is  not  very 
apparent  in  the  storage  modulus  plots  because  the 
system  frequency  difference  between  the  coated  and 
uncoated  specimens  is  very  small  near  the  transition 
region  and  the  modulus  values  calculated  from 
equations  (3.5)  and  (3.9)  are  prone  to  large  errors — 
as  much  as  50%  as  will  be  shown  in  a later  section. 

d)  The  measured  modulus  values  for  the  same 
coating  using  different  specimens  and  different 
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thicknesses  show  good  agreement.  Figures  (4.18)  and 
(4.19)  show  and  for  coating  #12325  as  functions 
of  temperature  at  two  different  frequencies 
calculated  from  two  torsion  pendulum  specimens  with 
different  coating  thicknesses.  Figures  (4.20)  and 

t »• 

(4.21)  show  E 2 and  for  the  same  coating  as 
functions  of  temperature  at  two  different  frequencies 
calculated  from  two  different  bending  specimens 
with  different  coating  thicknesses.  The  scatter  is 
greater  for  the  storage  modulus  data  because  of 
the  small  frequency  differences  used  in  its  computation. 
This  is  discussed  further  in  a later  section.  The 
consistency  with  which  the  low  temperature  values 
for  the  storage  modulus  for  a given  specimen  at  a 
specific  frequency  lie  above  or  below  the  trend 
line  is  because  smoothed  curve  values  were  used 
for  the  frequency  rather  than  individual  points. 

From  figures  (4.19)  and  (4.21)  it  can  be  seen 
that  the  maximum  values  of  G^  and  E^  for  coating 
#12325  are  0.74  x 10^  psi  and  1.34  * 10^  psi 
respectively. 

For  a material  with  stress  strain  relation 
o = M*e,  where  M*  = M'  + iM"  is  a complex  modulus, 
subjected  to  sinusoidal  deformation  c = cos  tut, 
it  can  be  shown  that  the  energy  dissipated  per  cycle 
is  given  by  [1] 


51 


AU  = TrM”e2 


(4.1) 


For  the  bending  specimen  configuration  with  anti- 
clastic  curvature  completely  constrained,  the 

2 

effective  Young's  modulus  is  given  by  E*/(l-v  ) i.e. 


ax  = (E*/(l-0)ex 


(4.2) 


so  that  the  energy  dissipated  is  given  by 


AU  = ffE" e2/ ( 1-v2) 


(4.3) 


where  is  now  the  amplitude  of  the  strain  in  the 

x direction  (figure  3.1).  From  Hooke's  Law: 

£ = (1/E)  (o  - vc  - vcr  ) 

x x y z 


e = (1/E)  (a  - va  - va  ) 
y y x z 

and 


e = (1/E)  (o  - vo  - vo  ) 
z / z x y 

With  no  anticlastic  curvature,  = 0 and  assuming 


oy  = 0,  we  get, 

ey  = - (v/(l-v) ) ex 


(4.4) 


From  a Mohr  strain  diagram,  it  can  be  seen  that  the 
amplitude  of  maximum  shear  stress  Ymax  is: 


Y = (l/(l-v))e. 
'max  0 


(4.5) 


Assuming  that  energy  dissipation  is  due  to  shear 


deformation  alone,  we  get 
AU  = ttg"y2 


max 

= TiG"e2/  (1-v)  2 


(4.6) 
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(4.7) 


From  equations  (4.3)  and  (4.6),  we  get 
e"  = ( ( 1+v ) / ( 1— v ) ) -G" 

M It 

Substituting  for  E and  G from  figures  (4.19)  and 
(4.21),  we  get 

1.34  = ( (l+v)/(l-v) ) 0.74 
or 

(l+v)/(l-v)  = 1.81 

which  corresponds  to  a value  of  v = 0.29.  The  value 

of  v for  coating  #16845  at  the  critical  temperature 

with  E"  = 1.30  x 106  psi  and  G"  = 0.53  x 106 
max  max 

psi  is  0.40.  Room  temperature  values  for  the  Poisson's 
ratio  for  glass  lie  between  0.2  and  0.25.  As  the 
temperature  increases,  the  value  of  v also  increases  [14]. 
At  the  firing  temperature  of  the  coating,  one  would 
expect  v to  have  reached  its  maximum  value  of  0.5; 
the  peak  in  the  loss  modulus  occurs  at  about  100  to 
200C  below  the  firing  temperature,  and  a value 
between  0.3  and  0.4  for  v at  this  temperature  is 
a reasonable  one. 

4.3.2  Distribution  of  Relaxation  Times 

The  method  of  reduced  variables  [35]  was  used  to 
analyze  the  temperature  and  frequency  dependence  of 
the  modulus  of  the  coating  material  in  the  transition 
region.  This  method  transforms  the  temperature  and 
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frequency  scales  so  that  the  observed  modulus  values 
can  be  brought  on  to  a single  curve  covering  a very 
wide  range  of  temperature  and  frequency.  It  is 
particularly  useful  if  the  available  frequency  range 
is  limited,  which  is  usually  true  for  simple  testing 
apparatus.  If  the  tests  can  be  made  at  various 
temperatures,  it  is  possible  by  this  method  to  predict 
the  effect  of  a much  wider  range  of  frequencies  than 
that  available  experimentally. 

The  procedure  was  as  follows.  First,  from  the 
graph  of  vs  temperature  at  3 different  frequencies 
(figure  (4.17)),  the  variation  of  with  the  logarithm 
of  frequency,  using  temperature  as  a parameter, 
was  obtained.  This  is  shown  in  figure  (4.22). 

II 

Figure  (4.23)  shows  this  variation  for  G^.  A 
reference  temperature  Tq  within  the  range  of 
experiments  was  chosen,  usually  about  the  middle  of 
the  transition  region.  At  each  temperature,  the 
average  frequency  factor  by  which  the  corresponding 
curve  (figure  (4.22))  had  to  be  shifted  to  match  the 
curve  corresponding  to  Tq  was  determined.  The 
logarithm  of  the  frequency  shift  factor  ntT  is  shown 
as  a function  of  T ^ in  figures  (4.24)  and  (4.25), 
for  the  bending  and  torsion  pendulum  cases 
respectively.  From  equation  (2.7),  this  curve  should 


The  points  fall  reasonably  close  to  a straight  line; 


so  the  best  straight  line  was  drawn  through  the 
points  by  inspection.  The  activation  energy,  Q,  of 
the  relaxation  process  was  determined  from  the  slope 
of  this  line  using  the  formula 

Q = 2.303  Rs  (4.8) 

where 

R = the  gas  constant 
= 1.987  cals/mole/°C 
and 

s = slope  of  the  line. 

The  E^  and  data  at  all  temperatures  was 

shifted  to  the  reference  temperature  by  using  the 

frequency  shift  parameter  aT.  (g")  was  found  from 

measurements  at  a given  temperature  T and  frequency 

f;  the  frequency  f was  multiplied  by  that  otT 

corresponding  to  T to  obtain  faT  = f^  which  is  the 

frequency  which  would  give  the  same  value  of  E2  (Gp 

at  the  reference  temperature  Tq.  Figures  (4.26) 

and  (4.27)  show  e"  and  g"  as  functions  of  frequency 

l 2 

for  coating  #12325  obtained  from  tests  on  specimens 
#16A  (Bending)  and  #2  (Torsion  Pendulum) . Some 
observations  on  the  dependence  of  E^  (G^)  on  frequency: 
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a)  The  points  in  figures  (4.26)  and  (4.27)  fall 
reasonably  close  to  a curve  that  peaks  symmetrically 
at  a certain  frequency.  The  curves  shown  have  been 
drawn  by  inspection. 

b)  The  width  of  the  peak  indicates  a distribution 
of  relaxation  times  rather  than  a single  relaxation 
time.  As  indicated  in  Chapter  (2.4),  the  width  at 
half  maximum  for  a single  relaxation  time  would  be 
equal  to  that  for  a Debye  peak,  i.e.  1.144  when 

E ^ (G^)  is  plotted  against  the  common  logarithm  of 
frequency. 

The  mean  relaxation  time  x , at  temperature  TQ 

was  deter:  : as  the  reciprocal  of  the  reduced 

circular  fr  . v , , at  which  e"  (G")  attains  a 

m 2 2 

maximum.  The  loss  modulus  w as  normalized  with 
respect  to  the  unrelaxed  modulus,  or  , and  the 
frequency  with  respect  to  the  frequency  f at  the 
peak  (corresponding  to  the  mean  relaxation  time  x ) . 
The  unrelaxed  modulus,  E^  or  G^,  was  chosen  as  the 
value  of  the  storage  modulus,  E'  or  G' , at 
approximately  150C  below  the  temperature  for  maximum 
loss  modulus.  This  is  approximately  the  temperature 
at  which  the  loss  modulus  begins  to  increase, 
indicating  the  beginning  of  the  relaxation  process. 
The  points  in  figures  (4.26)  and  (4.27)  were  then 
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replotted  on  the  curves  shown  in  figure  (2.5)  obtained 
from  a lognormal  assumption  for  the  distribution  of 
relaxation  times.  The  points  so  plotted  are  shown 
in  figures  (4.28)  and  (4.29).  By  inspection  of  the 
peak  height,  the  value  of  the  width  parameter,  B, 
was  obtained  as  6 for  the  bending  specimen  #16A, 
and  4.0  for  the  torsion  pendulum  specimen  #2.  For 
bending  specimen  #16A,  the  points  closely  follow  the 
curve  for  3=6.  For  torsion  specimen  #2,  the  points 
do  not  follow  the  curve  for  8 = 4.0  closely.  There 
is  more  scatter  in  the  torsion  pendulum  data. 

Figures  (4.30)  and  (4.31)  show  the  same  data  for  two 
other  bending  specimens  and  figures  (4.32)  and  (4.33) 
show  the  data  for  two  other  torsion  pendulum 
specimens.  The  bending  data,  in  general,  show  less 
scatter  than  the  torsion  pendulum  data.  For  bending 
specimens  #7  and  15A,  8 can  be  read  off  as  3 and  5 
respectively,  and  the  points  in  each  case  follow 
the  corresponding  curve  closely.  For  torsion 
pendulum  specimens  #1  and  2,  both  with  the  same 
coating,  #12325,  8 = 4 is  satisfactory.  The  points 
are  distributed  with  reasonable  scatter  around  the 
curve  for  6=4.  The  most  scatter  is  seen  in  the 
low  frequency  or  high  temperature  end.  For  torsion 
pendulum  specimen  #3,  the  peak  height  corresponds  to 
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8 = 5,  but  the  points  follow  a curve  approximately 
parallel  to  the  one  for  8=4.  From  the  evidence 
of  the  behavior  of  the  other  porcelain  enamels,  one 
can  disregard  the  possibility  that  the  enamel 
(#16845)  on  specimen  #3  does  not  obey  a lognormal 
relaxation  law.  The  discrepancy  observed  can  then 
be  ascribed  to  errors  in  measurement.  The  loss 
modulus  can  be  obtained  quite  accurately,  and  as 
shown  in  the  previous  section,  the  peak  value  of 
the  shear  loss  modulus,  obtained  from  the  torsion 
pendulum  experiments  compared  well  with  the  peak 
value  of  the  Young's  loss  modulus  obtained  from  the 
bending  experiments  for  the  same  coating,  with  a 
value  of  0.40  for  v at  the  transition  temperature. 
The  difficulty  is  in  getting  an  accurate  value  for 
the  unrelaxed  modulus  G^.  As  mentioned  earlier, 
(Chapter  4.3.1),  the  storage  modulus  G'  was  not 
calculated  for  the  coating  on  specimen  #3  because 
the  coating  was  so  thin  that  the  difference  in  the 
system  frequency  with  the  coated  and  uncoated 
specimens  was  extremely  small.  G^  was  assumed  to 
be  l/(2(l+v))  times  E^,  where  was  obtained  from 
a bending  test  on  the  same  coating  (specimen  #15A) . 
v was  assumed  to  be  0.25.  A change  in  the  choice  of 
v from  0.25  to  0.35  will  raise  each  point  in  figure 
(4.33)  by  as  much  as  8%.  The  points  need  be  raised 
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about  20%  to  match  with  the  curves  for  6=4.  There 
is  no  reason  to  believe  that  the  enamel  on 
specimen  #3  behaves  differently  from  other  enamels. 

It  is  believed  that  errors  in  measurements  caused 
the  difference  observed  in  the  case  of  specimen  #3, 
and  that  6 = 4 is  the  appropriate  value  for  that 
coating.  One  can  conclude  from  figures  (4.28)  to 
(4.33)  that  the  bending  measurements  are  more 
accurate  than  those  in  torsion. 

Tables  (4.7)  and  (4.8)  list  the  values  of  E~ 

z.  max 

(G2  max1 ' EU  (GU} ' fm  at  V T0 ' 6 and  Q for  the 
various  coatings  tested  in  bending  and  torsion.  For 

each  coating,  the  temperature  corresponding  to  a 

peak  frequency  f of  50  Hz  was  determined  using  the 

graph  of  log  otT  vs  T d for  that  coating.  The 

following  comments  can  be  made  about  the  values  shown 

in  Tables  (4.7)  and  (4.8). 

a)  For  a given  coating,  prepared  and  applied  under 
identical  conditions,  there  is  very  good  agreement 
in  the  values  of  the  loss  modulus,  the  temperature  Tq, 
the  width  parameter  8,  and  the  activation  energy  Q 
regardless  of  coating  thickness.  See  for  instance 
coating  S9919D  in  the  bending  specimens  and  coating 
#12325  in  the  torsion  specimens. 
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b)  For  the  same  coating  applied  under  different 
conditions  and  subjected  to  the  same  deformation, 
for  instance,  specimens  #16  and  16A,  and  specimens 
#15  and  15A  in  bending,  values  shown  in  the  Table  do 
not  agree.  The  coating  used  for  each  pair  of 
specimens  had  the  same  frit  number,  but  the  frit 
was  bought  in  different  batches.  Other  differences 
were  possibly  in  (1)  the  additions  to  the  frit-- 
normally  referred  to  as  mill  additions--used  to 
impart  certain  qualities  like  opacity,  ease  of  firing, 
etc.,  and  (2)  the  firing  process,  in  which  several 
factors  such  as  surface  preparation,  rate  of  heating, 
duration  of  hold  of  the  firing  temperature,  and  rate 
of  cooling  could  have  been  different.  The  coatings 
applied  at  the  University  of  Minnesota,  on  specimens 
#15A  and  16A,  tend  to  have  a higher  temperature  Tq,  for 
the  same  peak  frequency  f . As  has  been  discussed 
earlier,  it  is  suspected  that  the  high  value  for  the 
modulus  observed  for  the  coating  on  specimen  #16 

is  due  to  an  error  in  the  measurement  of  the  coating 
thickness . 

c)  For  coatings  prepared  and  applied  at  the  same 
place  but  subjected  to  different  states  of  stress 
(bending  specimen  #16A  and  torsion  pendulum  specimens 
#1  or  2) , there  is  a difference  in  the  value  of  the 
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width  parameter  6 measured.  From  figures  (4.28)  and 
(4.29)  it  can  be  seen  that  8 = 6 is  appropriate  for 
the  bending  case  while  8=4  fits  the  torsion  case. 
However,  the  values  of  E~  , g"  and  the 

temperature  for  the  same  peak  frequency  agree 
very  well.  The  width  of  the  peak  for  the  loss 
coefficient  vs.  temperature  curve  was  also  narrower 
for  the  torsion  case  (Tables  (4.2)  and  (4.3)).  Since 
system  measurements  depend  on  integrated  properties 
through  the  whole  volume  being  tested,  the  apparent 
dependence  on  temperature  is  reduced  if  different 
parts  of  the  specimen  are  at  different  temperatures. 

Since  the  temperature  distribution  in  the  torsion 
specimen  is  less  uniform  than  that  in  the  bending 
specimen,  any  errors  in  the  temperature  measurement 
should  cause  a wider  peak  in  the  torsion  case  than 
in  bending.  Hence,  we  can  safely  say  that  the  difference 
in  8 between  the  two  cases  does  not  arise  from 
erroneous  temperature  measurements.  The  lower  8 
in  torsion  is  accompanied  by  a lower  measured 
activation  energy  Q.  The  difference  in  8 between 
torsion  and  bending  cases  is  also  apparent  for  the 
coating  #16845. 

There  appear  to  be  two  possible  reasons  for  the 
difference  in  8.  1)  The  state  of  stress  to  which  the 
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coating  is  subjected,  sinusoidal  shear  in  one  case 
(superimposed  upon  a very  small  static  tension  due 
to  the  weight  of  the  pendulum)  and  constrained 
uniaxial  tension  in  the  other.  2)  The  composition 
of  the  base  metal  in  the  two  cases.  The  torsion 
pendulum  specimens  were  made  of  multimet  alloy 
due  to  lack  of  immediate  availability  of  Hastelloy 
X stock  of  the  required  size. 

In  order  to  investigate  this,  a bending 
specimen  (#M1)  was  made  of  Multimet  alloy  and 
coated  with  frit  #12325.  Tests  on  this  specimen 
gave  the  following  results. 

Ey  - 9 x 106  psi 

E"  = 1.7  x 106  psi 

2 max 

Tq  (for  fm  = 50  Hz)  = 589C 

Q = 220  kcals/mole 

and 

6 = 4.5 

The  modulus  values  agree  well  with  the  previous 
measurements  on  this  coating.  The  temperatures 
measured  for  this  specimen  were  in  general  about  15C 
lower  than  those  for  specimen  #16A  or  for  the  torsion 
specimens.  However,  what  is  interesting  is  that  the 
values  for  activation  energy  Q and  the  width 
parameter  0 agree  better  with  the  values  obtained 


the  case  of  specimen  Ml. 

d)  The  value  of  8 does  not  have  any  observable 
relation  with  the  firing  temperature  of  the  coating. 

e)  An  increase  in  8 is  associated  with  an  increase 
in  Q,  i.e.  the  wider  the  distribution,  the  greater 
the  activation  energy. 

f)  The  activation  energies  measured  are  very  high. 

The  activation  energy  for  a diffusion  mechanism  is 
typically  30  kcals/mole.  The  activation  energy  for 
coating  S9919D  and  #16845  are  within  the  range  of 
activation  energies  for  viscous  flow  in  glass  forming 
silicon  dioxide  [12] . But  the  activation  energy  of 
about  300  kcal/mole  for  coating  #12325  is  considerably 
higher . 

Figures  (4.34)  and  (4.35)  show  the  values  of 
and  e"  calculated  from  the  actual  measured  data 
for  specimen  #2  in  torsion  and  specimen  #16A  in 
bending  compared  with  the  curve  calculated  from  the 
lognormal  distribution  assumption  using  the  parameters 
determined  as  shown  above.  The  agreement  is  well 
within  experimental  error,  particularly  in  the 
temperature  region  where  the  loss  modulus  peaks.  (An 
error  analysis  is  made  in  the  next  section.) 

In  determining  the  frequency  shift  factor  aT  as 
a function  of  temperature,  the  loss  modulus  rather 


than  the  storage  modulus  was  used  because  the  error 


64 


involved  in  the  calculation  of  the  loss  modulus  is 
smaller  and  there  is  less  scatter  in  the  data.  This 
can  be  seen  in  figures  (4.18),  (4.19),  (4.20)  and 

(4.21).  Once  aT  is  known  as  a function  of  temperature, 
it  can  be  used  to  shift  the  storage  modulus  data  also 
to  a reference  temperature.  The  points  taken  from 
a smoothed  curve  of  e!,  vs.  temperature  and  vs. 
temperature  were  reduced  to  temperature  Tq  and 
normalized  and  are  shown  in  figures  (4.36)  and  (4.37) 
plotted  with  the  normalized  curves  of  figure  (2.6) 
obtained  from  a lognormal  assumption  for  the 
distribution  of  relaxation  times.  Figures  (4.38) 
and  (4.39)  show  the  actual  measured  values  of 
and  E ^ compared  with  the  curve  of  G^  or  E^  vs. 
temperature  calculated  from  the  lognormal  assumption 
using  the  parameters  determined  above.  It  can  be 
seen  that  the  agreement  is  again  within  experimental 
error. 

4.4  Error  Analysis 

The  storage  modulus  E^  and  the  loss  modulus  E^ 
were  calculated  from  measured  data  using  the 
equations 


- <VC)  <£c/fb  - 

(4.9) 

and 

E2  - <V°>  * fc/fb 

(4.10) 
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where 


E^  - Young's  modulus  of  base  metal 
c - 3t  + 3t2  + t3 
t = thickness  ratio 

f = system  frequency  with  coated  specimen 
fb  = system  frequency  with  bare  specimen 
and 

n = system  loss  coefficient  due  to  the 
coating  alone. 

This  method  of  calculation  introduces  large 
errors  in  the  value  of  even  for  small  errors  in 
the  measured  frequency  values.  It  can  be  seen 
that  the  calculation  of  the  storage  modulus  depends 
critically  upon  the  small  differences  in  the  system 
frequency  with  the  coated  and  uncoated  specimens. 

Let  A ( ) represent  the  relative  error  in  the 

quantity  in  parentheses.  Assuming 
A (f ) = ±ef 

A(fc/fb)  = ±2ef 
and 

4l£c/fb>  - t4ef 
Therefore 

A(f2/f2  - 1)  = (±4ef)/(f2/f2  - 1) 
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In  the  transition  region,  f /f^  is  about  1.04. 
therefore, 

A(fc/fb  " 1)  = ±100ef  (4.11) 

The  thickness  of  the  coating  cannot  be  controlled 
accurately  while  applying  it.  It  has  therefore  to 
be  measured  after  application.  It  was  obtained  as 
the  difference  in  specimen  thickness  (diameters) 
before  and  after  application  of  the  coating.  Let  d.  * 
be  the  coated  specimen  thickness 

d^  is  the  bare  specimen  thickness 
and 

d2  = d^*  - d^  is  the  coating  thickness. 

Suppose  the  relative  error  in  measuring  d is  e^ . 
i.e.  A (d)  = ±ed 
Therefore , 

Afd^  = ±ed 
and 

A^*)  = ±ed 

Then, 

A(d2)  = A(dd*  - dd)  = +2ed  dd/d2 
= ±40ed  with  ~ ^0. 

Therefore , 

A(d2/dd)  = A(t)  = ±40ed. 

2 3 

Now,  c = 3t  + 3t  + t , and  because  t is  very  small 
compared  to  1,  we  get 
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A (c)  = ±400^  = ±ec  (4.12) 

Now, 

i(Ej)  = A(f*/f£  “ 1}  + A(c) 

= ±100  e-  ± 40  e,  (4.13) 

f a 

It  is  estimated  that  in  the  worst  case,  the 
error  in  the  frequency  measurements  can  be  controlled 
to  0.3%  and  the  error  in  measuring  the  specimen 
thickness  to  0.05%  (Appendix  (6)).  A conservative 
estimate,  which  allows  for  other  inaccuracies  that 
might  occur  during  the  reduction  process,  as  when 
matching  bare  specimen  curves,  is  that  the  elastic 
modulus  of  the  coating  material  can  be  determined 
from  equation  (4.9)  using  measured  frequency  and 
thickness  values  with  an  error  of  less  than  50%. 

The  least  accurate  values  are  at  temperatures  above 
the  transition  temperature  region  where  the 
difference  between  the  frequencies  with  the  coated 
and  uncoated  specimens  is  very  small  and  measurement 
errors  are  magnified  considerably. 

Now  from  equation  (4.10), 

AfE^)  = ±4ef  ± 40ed  ± A(n). 

The  error  in  measuring  n is  relatively  low.  By 
subtracting  out  the  system  loss  coefficient  with  the 
specimen  uncoated,  not  only  the  base  material  damping. 
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but  also  other  extraneous  losses  like  air  friction 
and  support  effects  are  eliminated.  Also,  the 
measured  system  loss  coefficient  with  the  bare 

specimen  in  the  temperature  range  of  interest  is 

- 3 

very  low,  about  10  , (except  for  the  high  frequency 

torsion  specimen) , while  the  loss  coefficient  with 

the  coating  in  the  transition  region  in  most  cases 
-2 

exceeds  10  . The  loss  modulus,  therefore,  can  be 

calculated  from  measured  data  using  equation  (4.10) 


with  a maximum  error  of  ±5%. 


Chapter  5.  Conclusions  and  Discussion 


A test  program  designed  to  study  the  damping  to 
be  had  from  porcelain  enamel  coatings  at  temperatures 
where  they  soften  was  described  in  the  preceding 
chapters.  Three  different  pieces  of  apparatus  were 
built;  the  first  measured  the  vibration  decay 
of  a coated  solid  cylindrical  specimen  in  torsion  in 
the  frequency  range  1 Hz  to  15  Hz.  The  second 
measured  the  decay  of  a coated  flat  plate  specimen 
in  bending  in  the  frequency  range  11  Hz  to  170  Hz, 
and  the  third  measured  the  decay  of  a coated  hollow 
cylindrical  specimen  in  torsion  in  the  frequency 
range  165  Hz  to  600  Hz.  Several  different  enamel 
frits  were  used  in  the  tests;  the  lowest  effective 
damping  temperature  was  400C  and  the  highest  was 
750C.  The  coating  thicknesses  ranged  from  0.001  to 
0.006  inches.  From  the  test  results  presented  in 
Chapter  4,  the  following  broad  conclusions  can  be 
drawn  about  the  equipment  used  and  about  the  porcelain 
enamels . 

1)  The  torsion  pendulum  and  the  cantilever  beam 
apparatus  can  be  used  satisfactorily  to  obtain  the 
complex  modulus  of  porcelain  enamels  as  a function 
of  temperature. 
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2)  The  hollow  cylindrical  torsion  specimen  apparatus 
is  not  a satisfactory  piece  of  equipment  because  of 
the  expansion  joints  in  the  specimen  assembly  which 
caused  high,  non-linear  damping.  The  loss  coefficient 
data  obtained  from  this  equipment  in  the  frequency 
range  165  to  600  Hz  was  not  used  to  calculate 
the  complex  modulus;  but  it  does  show  the  same 
general  trend  as  the  data  at  the  lower  frequencies. 

3)  The  damping  tests  on  some  specimens  were 
repeated.  The  results  obtained  were  the  same  as 
for  previous  tests.  This  indicates  that  the  testing 
procedure  does  not  damage  the  coating,  that  the 
heating-cooling  cycle  does  not  change  the  mechanical 
properties  of  the  enamel,  and  that  the  relaxation 
process  is  reversible. 

4)  The  interface  between  the  enamel  and  the  metal 
is  not  important  as  there  is  very  good  agreement 
among  tests  using  different  enamel  thicknesses. 

5)  The  storage  modulus  for  the  various  enamels  is 
approximately  the  same.  At  room  temperature,  the 
storage  Young's  modulus  is  about  10  x 10^  psi,  and 
the  shear  modulus  measured  for  one  frit  (#12325)  is 
about  5 x io^  psi.  These  agree  well  with  values 
for  glass  reported  elsewhere. 

6)  The  loss  modulus  shows  a symmetric  peak  when 
plotted  against  either  temperature  or  frequency. 
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The  maximum  value  is  1.4  x 10^  psi  for  the  loss 
Young's  modulus  ar.d  0.7  * 10^  psi  for  the  loss  shear 

• odulus . 

The  widths  of  the  peaks  in  the  loss  modulus  vs. 
ature  curves  indicate  multiple  relaxations. 

Ti  lognormal  distribution  of  relaxation  times  was 
u 1 to  describe  the  relaxation  process.  The  scatter 
i..  : he  experimental  data  is  high,  because  the 
porcelain  enamels  can  be  tested  only  as  thin 
coatings,  and  a lognormal  distribution  function  for 
the  relaxation  process  is  an  adequate  one.  Three 
parameters  suffice  to  describe  the  behavior  of  the 

" II 

porcelain  enamels,  namely#  E^  ax  or  G 2 miy>  S and 
t . Since  the  relaxed  modulus,  E0  or  can  be 

in  k iv. 

assumed  to  be  zero,  the  unrelaxed  modulus,  E^  or 
Gl,,  8 and  t can  be  used  to  completely  describe 
the  relaxation  process. 

8)  The  temperature  bandwidth  varies  between  50C 
and  70C  for  all  enamels. 

9)  The  frequency  bandwidth,  6#  is  different  for 
different  enamels  ar.d  has  no  observable  relationship 
to  the  softening  temperature  of  the  porcelain 
enamel . 

10)  The  composition  of  the  enamel  significantly 
affects  the  softening  temperature  and  the  width  8, 
but  not  the  modulus.  The  0's  measured  for  the 
same  coating  on  different  base  metals  are  also 


different.  It  is  possible  that  during  the  initial 
firing  process,  there  is  some  diffusion  of  ions 
between  the  enamel  and  the  base  metal.  The  coating 
being  very  thin,  this  would  significantly  alter  its 
composition;  that  is,  the  composition  of  the  coated 
enamel  xs  different  from  that  of  the  uncoated  frit. 
The  final  composition  of  the  coating  would  then 
depend  on  the  composition  of  the  base  roe  cal  also. 
Such  an  interaction  must  be  complete  during  the 
firing  process,  because  later  heating-cooling  cycle-: 
do  not  change  the  measured  properties  of  the 
enamel.  Also,  the  interface  between  the  enamel 
the  metal  surface  is  not  important. 

11)  The  relaxation  process  follows  an  Arrhenius 
equation  for  several  decades  of  frequency- -at  ic 
7 decades  in  bending  (figure  (4.24)),  and  at  least 
5 decades  in  torsion  (figure  (4.25)).  This  does  not 
imply  that  the  enamels  are  Arrhenius  over  a smaller 
frequency  range  when  subjected  to  torsional  strain. 
The  measured  range  of  Arrhenius  behavior  is  smaller 
in  torsion  only  because  i)  the  torsion  measurements 
were  made  over  a smaller  frequency  range;  they 
covered  a factor  of  11.6  compared  with  a factor  of 
15.5  in  bending,  and  ii)  the  activation  energy 
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measured  in  bending  is  greater  than  that  in  torsion. 
It  has  been  pointed  out  earlier  that  this  did  not 
depend  on  the  state  of  stress,  but  was  due  to  the 
base  metals  in  the  two  cases  being  different,  thus 
altering  the  final  composition  of  the  coating. 

12)  The  activation  energies  measured  are  very  high, 
of  the  order  of  200  kcals/mole.  The  typical 
activation  energy  for  alkali  ion  diffusion  in  glass 
is  20  - 30  kcals/mole;  that  for  viscous  flow  is 
known  to  be  higher,  about  100  kcals/mole.  Even  the 
energy  of  the  silicon-oxygen  bond  is  only  109 
kcals/mole.  Extremely  high  values  for  activation 
energy,  between  120  and  190  kcals/mole,  have  also 
been  measured  for  stress  relaxation  in  sodium 
phosphate  glasses  in  the  transition  temperature 
region  [20].  Very  high  activation  energies,  up  to 
242  kcals/mole,  have  been  reported  for  viscous 
flow  in  alkali  borate  glasses  [45] . Polymers  in  the 
transition  region  also  exhibit  such  high  activation 
energies  [46]  . No  satisfactory  explanation  is 
available  for  these  high  activation  energies.  A 
value  of  200  kcals/mole  cannot  be  associated  with 
any  known  physical  mechanism.  One  can  speculate  that 
the  network  structure  of  the  enamel  is  such  as  to 
impose  severe  constraints  on  the  atomic  movements 
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which  cause  the  relaxation  in  the  material,  thus 
increasing  the  activation  energy  of  the  relaxation 
process.  One  can  also  speculate  that  the  atomic 
movement,  when  it  does  take  place,  cannot  do  so 
without  disturbing  other  neighboring  atoms,  thus 
initiating  a sequence  of  motions.  It  is  apparent 
that  such  a cooperative  process  would  be  extremely 
complex  and  cannot  be  described  in  detail.  It 
is  also  doubtful  that  one  can  assign  a single 
activation  energy  to  such  a complicated  process. 

This  is  the  first  work  that  investigates  the 
damping  properties  of  porcelain  enamels  in  detail. 
It  has  left  unanswered  some  questions  regarding 
1)  the  effect  of  the  coating  composition  and  2) 
the  anomalously  high  activation  energies,  which 
provides  interesting  areas  for  future  work.  It  has 
demonstrated  convincingly  that  porcelain  enamels 
can  be  used  as  efficient  damping  materials,  and 
that  there  is  a need  for  better  and  more  accurate 
testing  equipment.  It  would  be  worthwhile  to 
develop  porcelain  enamels  with  lower  softening 
temperatures.  Porcelain  enamels  with  softening 
temperatures  up  to  about  1200C  are  being  used;  but 
there  are  hardly  any  with  softening  temperatures 
lower  than  400C.  It  would  be  useful  to  extend  the 
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temperature  range  over  which  porcelain  enamels,  as  a 
class,  can  be  used  efficiently  for  damping  purposes. 


Appendix  1.  Porcelain  Enamels 

Porcelain  enamels  are  thin,  glassy,  inorganic 
coatings  bonded  to  a metal  surface  by  fusion.  They 
are  usually  laid  in  thicknesses  varying  from 
0.001  inch  to  0.01  inch.  They  are  primarily  oxides 
and  are  highly  resistant  to  oxidation  and  corrosion. 
They  are  smooth  and  hard  and  very  resistant  to 
thermal  stresses  and  shock.  Porcelain  enamels  are 
extensively  used  to  provide  protection  for  metal 
parts  exposed  to  corrosive  and  erosive  environments. 
They  behave  anelastically  in  their  softening 
temperature  range  and  hence  are  good  candidate 
materials  for  damping  at  high  temperatures. 

Composition  and  Method  of  Coating  [6] 

Porcelain  enamels  consist  of  raw  materials  such 
as  quartz,  feldspar,  clay,  fluorspar,  borax,  soda 
ash,  kryolite,  nitrate  of  soda,  zirconium  oxide, 
titanium  oxide,  cobalt  oxide,  nickel  oxide,  and 
manganese  oxide  which  are  fused  together  to  form 
the  frit  in  granular  form.  The  frit  is  finely  powdered 
along  with  some  additional  chemicals  to  impart  certain 
desirable  properties  to  the  enamel.  The  chemicals 
added  include  floating  agents,  opacifiers,  coloring 
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agents  and  electrolytes  to  deflocculate  the  clay  and 
properly  suspend  the  enamel.  The  powder  is  then 
mixed  thoroughly  with  water  as  the  vehicle  to  form 
a slurry  or  slip,  and  is  ready  for  application. 

The  metal  surface  to  be  coated  is  prepared  by 
chemical  methods  (like  acid  pickling)  or  by  mechanical 
methods  (like  sand  blasting)  to  produce  a roughened 
surface  free  from  any  metallic  oxides  or  organic 
material.  The  slip  is  then  applied  to  the  clean 
surface  by  dipping,  brushing  or  spraying,  and  dried 
in  air.  The  dried  coating  is  fired  for  two  to 
fifteen  minutes  in  an  oxidizing  atmosphere  at  the 
temperature  required  to  re-fuse  the  glassy 
constituents,  until  the  coating  material  bonds  with 
the  metal  to  form  a homogeneous  vitreous  coating  [36], 
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Appendix  2.  Relaxation  in  Porcelain 


A2 . 1 Mechanical  Models 

The  most  general  form  of  a linear  constitutive 
equation  to  describe  an  anelastic  solid  is  [40] 


(A2.1) 

where  the  a. 's  and  b.'s  are  constants,  a is  the 

l i 

stress  and  e is  the  strain.  The  behavior  of  an 
anelastic  material  is  often  described  by  mechanical 
models.  A mechanical  model  is  appropriate  to  a 
particular  material  if  the  relations  between  force, 
displacement  and  time  for  the  model  are  the  same  as 
those  for  stress,  strain  and  time  in  the  material. 
There  are  two  basic  elements  in  any  mechanical 
model — the  ideal  spring  and  the  Newtonian  dashpot. 

The  ideal  spring  alone  represents  an  elastic  body 
for  which  only  Sq  and  b^  in  equation  (A2.1)  are 
different  from  zero.  In  the  dashpot  alone,  the  force 


, . d d , . 

(ao  ai  dt  a2  771  * * ,)0 

at 


or 


“ ,bo  + bi  3t  + b2  TV  + 

dt 


a1 

( E a.  -SU-)ct  = ( E b.  — r)e 
i=0  dt  i=0  dt 


• . ) e 


is  proportional  to  the  velocity,  or,  only  ag  and  in 
equation  (A2.1)  are  different  from  zero.  An 
anelastic  solid  is  represented  by  a combination  of 
springs  and  dashpots--the  springs  providing  the 
elasticity  and  the  dashpots  the  internal  friction. 

The  simplest  models  that  can  be  constructed 
from  the  two  basic  elements  are  the  Voight  and  the 
Maxwell  models.  The  Voight  model  has  a spring  and 
dashpot  connected  in  parallel  while  the  Maxwell 
model  has  the  elements  connected  in  series.  These 
two  are  rather  too  simple  to  represent  any  real 
material.  The  simplest  model  which  is  representative 
of  a real  material  is  what  is  called  the  standard 
three  parameter  modal,  consisting  of  a Voight  unit 
in  series  with  a spring.  It  can  be  shown  [41]  that 
for  a sinusoidal  deformation,  the  modulus  M*  = M'  + iM" , 
where  M*  is  the  Young's  modulus  for  longitudinal 
deformation  or  the  shear  modulus  for  shear  deformation, 
is  given  by 

M'  (to)  = Md  + 6M  (<o2t?)/(1  + to2T2)  (A2.2) 

and 

M"  (to)  = 6M  toT£/(l  + u)2T2)  (A2.3) 

where 

My  = unrelaxed  modulus 

M = relaxed  modulus 
K 


80 


6M  = M - M = magnitude  of  the  relaxation 

U t\ 

w = circular  frequency 

x£  = relaxation  time  at  constant  strain 
= y6M/M~ 

and 

M = viscosity  of  the  dashpot. 

M'  decreases  from  M to  M as  the  frequency  decreases. 
The  function  M"  reaches  its  maximum  value  at  a 
frequency  such  that  cox£  = 1.  When  plotted  against 
log  (ux£) , the  curve  has  a width  at  half  maximum 
of  1.144.  This  is  called  a Debye  peak. 

A2 . 2 Lognormal  Distribution  of  Relaxation  Times 

For  an  anetnstic  material  with  a distribution 

function  ( iln  x)  for  the  logarithm  of  the  relaxation 

times  x,  the  storage  modulus  M'  and  the  loss  modulus 

M"  are,  as  shown  in  equations  (2.1)  and  (2.2) 

M'  (w)  + M_,  + 6M/°°  ( 4>  ( £n  x)  w2x2)  / (l+w2x  2)  d£n  x 
K 

(A2.4) 

and 

oo 

M"  (to)  = 6M/  (<^(in  x)  wx)/(l+w2x2)djln  x (A2.5) 

— OO 

where 

00 

/ <{>(in  x)d£n  x = 1 (A2.6) 

— 00 
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Assuming  a normalized  Gaussian  distribution  in  An  t, 

and  introducing  a new  variable  z = in(i/T  ) , where 

t is  the  mean  value  of  x,  we  get 
m 3 


<Hz)  = ( 1/  ( B /tT)  ) exp.  (-z2/g2)  (A2.7) 

where  8 is  the  value  of  |z|  for  which  4>(z)  falls  to 
1/e  of  its  maximum  at  z = 0 . Then, 

oo 

M"  (w)  = 6M  ; (<p  (z)wx)/(l  + U2x2)dz  (A2.8) 

— OO 

Introducing  the  variable  x = £n(ojxm), 

00  ? 

M"(x)/6M  = (l/(2/r F)  ) / exp  (-u  ) sech  (x+8u)du 

— CO 

= f (x  , 8 ) (A2.9) 

Thus , 

M"  (x  ) /M"  ( 0 ) = f (x1  ,8)/f  (0,6)  = F (x  ' , 8 ) (A2.10) 

where  x'  = log  (ux^)  = x/2.303.  Also, 

M" (0)/6M  = f (0, 8)  (A2. 11) 

Reference  [40]  tabulates  the  function  F for  different 
values  of  x'  and  8,  and  2f(0,B)  for  different  B's. 

A portion  of  these  tabulations  are  shown  in  Tables 
(2.2)  and  (2.3).  Figure  (2.3)  shows  plots  of  F vs  x' 
for  various  values  of  8.  The  case  8=0  corresponds 
to  the  case  of  a material  with  a single  relaxation 
time.  For  all  8,  the  function  F(x',B)  reaches  its 
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maximum  value  of  1 at  = 1 or  x'  = 0.  Figure 

(2.4)  shows  2 f (0r@)  as  a function  of  3.  The  function 
f(x',B)  was  obtained  from  F(x’,3)  and  f(0,8)  and  is 
tabulated  in  Table  (2.4).  Figure  (2.5)  shows  plots 
of  f(x',3)  vs.  x'  for  various  8's.  It  can  be  seen 
from  figures  (2.4)  and  (2.5)  that  the  peak  value 
of  the  loss  modulus  decreases  with  increasing  3. 

The  storage  modulus  M'  is  given  by  equation 
(A2.4).  For  porcelain  enamels,  we  shall  assume  that 
the  relaxed  modulus  MR  is  very  small,  and  that 
6M  = My.  Therefore, 

00 

M 1 (to)  = M / (<j>Un  t ) uj2t2)  / (1  + cu 2 t 2 ) din  x (A2.12) 


with  the  distribution  given  by  equation  (A2.7),  and 
for  z,  x,  and  x'  defined  as  above,  we  get 


(M'/My  - 0.5)  = £1  (x' ,3) 


(A2.13] 


where 


f^(x',8)  = (1/(2/tt))/  exp  (-u  ) tanh(x+3u)du 


(A2.14) 


It  can  be  shown  that  for  all  8, 


f1(0,g)  = 0 


(A2.15) 


M ' /MtI  = 0 . 5 at  x = 0 or  WT  = 1. 
u m 


Table  (2.5)  shows  f^(x',6)  for  different  values  of 
x'  and  3.  The  integral  defined  in  equation  (A2.14) 


Ml1  1ST 


was  evaluated  numerically.  Only  positive  values  of 
x'  are  shown  as  the  function  f-(x',8)  is  odd 
about  x'  =0.  Figure  (2.6)  shows  plots  of  M'/My 
as  a function  of  x'  for  various  values  of  8. 


This  apparatus  was  used  in  the  frequency  range 
1 Hz  to  15  Hz.  The  solid  cylindrical  test  specimen 
shown  in  figure  (A3.1)  was  made  of  multimet  alloy, 
(Table  (4.9)),  15  inches  long  and  0.1  inch  in 
diameter  with  flat  plates  welded  to  the  ends  to  permit 
clamping.  The  middle  4 inches  of  the  specimen  was 
coated  with  the  porcelain  enamel.  The  apparatus  used 
for  the  tests  is  shown  in  figure  (A3. 2).  The  specimen 
was  rigidly  clamped  at  the  top  end  to  a massive  block. 
The  3/16  inch  diameter,  1/2  inch  projection  at  the 
lower  end  was  fastened  to  a cylindrical  steel 
adapter,  the  lower  end  of  which  had  a 1/8  inch 
projection  that  was  guided  through  a simple  ball 
bearing.  A steel  driving  arm  carrying  an  aluminum 
capacitance  plate  was  fastened  to  the  lower  plate. 

The  system  could  be  set  into  free  vibration  at  9.0 
Hz  by  plucking  the  arm.  The  frequency  of  vibration 
was  reduced  to  1.3  Hz  at  room  temperature  by  clamping 
the  mass  shown  in  figure  (A3. 3)  to  the  lower  plate. 

A frequency  of  15.0  Hz  at  room  temperature  was 
obtained  by  using  an  aluminum  driving  arm  and 
capacitance  plate.  The  capacitance  plate  moved 
between  two  pairs  of  differentially  connected  plates, 
and  the  change  in  capacitance  caused  by  the  change 
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in  area  covered  by  the  plate  as  the  specimen  vibrated 
was  used  to  generate  a voltage  proportional  to  the 
vibration  amplitude. 

A cylindrical  electric  furnace  with  three  zones 
with  separate  resistance  controls  was  used  to  heat 
the  specimen  and  hold  it  at  the  desired  temperature. 
The  furnace  had  a bore  of  2.5  inches  and  was  12 
inches  long,  so  that  1.5  inches  of  the  specimen 
projected  outside  the  furnace  at  each  end.  During 
the  tests,  the  furnace  chamber  was  closed  by  insulat- 
ing plates.  A chromel-alumel  thermocouple,  its  lead 
wires  enclosed  in  refractory  tubes,  inserted  into  the 
furnace  chamber  through  a hole  in  the  furnace  wall 
served  to  control  the  furnace  temperature  and  was  also 
the  indicating  thermocouple  during  the  tests.  It 
could  be  moved  radially  to  be  within  0.03  inches  of 
the  specimen  surface.  A temperature  survey  was 
conducted  with  8 chromel-alumel  thermocouples  fastened 
at  various  points  along  the  specimen  length,  as  shown 
in  figure  (A3. 4).  The  zone  resistors  were  adjusted 
so  that  the  temperature  across  the  test  section  of 
the  specimen  was  uniform  within  20F  (11C) . The 
temperature  indicated  by  the  control  thermocouple  was 
found  to  be  20F  (11C)  less  than  the  actual  temperature 
of  the  specimen. 
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A block  diagram  of  the  instrumentation  used  is 
shown  in  figure  (A3. 5).  The  oscillatory  change  in 
capacitance  was  converted  to  an  analogous  voltage 
signal  using  a suitable  electronic  transducer.  A 
band  pass  filter  was  used  to  filter  out  any  stray 
signals — like  high  frequency  electronic  noise — and 
give  a clean  signal  which  was  monitored  on  an 
oscilloscope.  A level  recorder  which  gave  the 
envelope  of  the  decay  curve  to  a logarithmic  scale 
was  used  at  the  higher  frequencies  of  9 and  15  Hz. 

A typical  curve  is  shown  in  figure  (A3. 6).  The 
slope  of  this  curve  is  directly  proportional  to  the 
loss  coefficient.  For  tests  made  at  1.2  Hz,  the 
voltage  signal  was  recorded  directly  using  a pen 
recorder  and  the  loss  coefficient  determined  from 
a semi-logarithmic  replot  of  the  data.  The 
frequency  counter  could  be  used  for  a 1,  10  or  100 
second  measurement  and  by  taking  an  average  of 
several  measurements  the  frequency  could  be 
obtained  to  within  ±0.1%. 


Natural  Frequency 


The  strain  energy  for  data  reduction  was  calcu- 
lated from  an  assumed  stiffness  of  the  specimen. 

This  stiffness  was  verified  by  comparing  the 
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measured  and  the  computed  natural  frequencies  of  the 
system. 

For  a uniform,  solid,  cylindrical  rod,  the 
torsional  stiffness,  K,  is 


K = ?fi G 

where 

d = diameter  of  rod 
l - length  of  rod 


G = shear  modulus  of  the  material  of  the  rod. 
For  d = 0.1  inch,  H = 15  inches  and  G for  mul timet 

g 

alloy  at  room  temperature  = 11.6  * 10  psi, 

K = 7.6  inch  pounds/radian. 

The  frequency  of  a torsion  pendulum  is  given  by 


f = JZ 

1 2tt  / J 


0.44 


where  J = mass  moment  of  inertia  of  the  pendulum  bob. 
The  computed  values  of  J,  the  corresponding 
frequencies  and  the  measured  room  temperature 
frequencies  were  as  follows. 


with 

J 2 
in  lb.  sec 

f Hz 

computed 

f Hz 

measured 

Aluminum  driving  arm 
and  capacitance  plate 

0.77X10-3 

15.8 

15.5 

Steel  driving  arm  and 
capacitance  plate 

2.4xl0-3 

9.0 

8.99 

Steel  driving  arm  and 
capacitance  plate  and 
steel  weight  of 
figure  (A3. 3) 

0.114 

1.303 

1.305 

Since  the  measured  natural  frequencies  for  all 


three  configurations  agreed  well  with  the  computed 
frequencies  using  a stiffness  value  of  7.6  inch 
pounds/radian,  this  value  of  stiffness  was  used  for 
strain  energy  calculations  during  data  reduction. 


Appendix  4.  Cantilever  Beam  Apparatus 

The  original  bending  test  apparatus  was  designed 
by  B.  J.  Lazan  and  was  modified  by  Lee  [8]  to  operate 
at  32  Hz.  Lee's  apparatus  is  shown  in  figure  (A4.1). 
It  was  later  modified  as  shown  in  figure  (A4.2)  to 
operate  at  47  Hz.  The  flat  test  beam  made  of 
Hastelloy  X (Table  (4.9))  was  firmly  fastened  at 
one  end  to  a support  which  was  tied  down  to  a massive 
base.  An  arm  made  of  Hastelloy  X,  extending  down- 
wards outside  the  furnace  was  fastened  to  the  free 
end  of  the  specimen  for  measurement  and  control 
outside  of  the  furnace.  The  arm  could  be  displaced 
manually  and  a hook  at  its  lower  end  held  by  a pawl 
carried  at  the  end  of  a micrometer  screw  which  was 
used  to  adjust  the  initial  displacement.  The  system 
was  set  into  free  vibration  at  its  fundamental  natural 
frequency  by  releasing  the  arm.  The  natural 
freqeuncy  of  vibration  of  the  system  could  be 
decreased  by  clamping  steel  weights  to  the  lower  end 
of  the  arm.  The  lowest  freqeuncy  so  obtained  using 
reasonably  sized  weights  was  11  Hz.  A piezoelectric 
accelerometer  fastened  to  the  arm  outside  the  furnace 
was  used  to  pick  up  the  vibration  signal. 

The  test  specimen,  shown  in  figure  (A4.3)  had 
the  same  geometry  as  the  one  used  in  the  original 
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indicated  the  temperature  of  the  specimen  during  the 
tests.  For  each  test,  the  temperature  should  be 
uniform  throughout  the  specimen  and  constant  in  time. 
A soak  time  test  was  performed  to  determine  the  time 
taken  to  reach  steady  temperature  and  a temperature 
survey  was  conducted  to  measure  the  temperature 
variation  along  the  specimen.  The  specimen  used 
for  this  survey  had  two  chromel-alumel  thermocouples 
fastened  to  each  face;  other  chromel-alumel  thermo- 
couples were  located  as  shown  in  figure  (A4.4). 
Thermocouple  #7  was  used  to  monitor  the  temperature 
at  the  point  where  the  accelerometer  was  fastened. 
Since  the  maximum  operating  temperature  of  the 
accelerometer  used  in  some  tests  was  350F  (~180C) , 
a small  blower  was  used  to  keep  its  temperature 
within  safe  limits.  It  was  found  that  it  took  2 1/2 
hours  to  reach  temperature  equilibrium  for  the  first 
relatively  large  temperature  change,  and  45  minutes 
for  each  subsequent  smaller  temperature  increment. 

The  temperature  variation  between  the  upper  and  lower 
sides  of  the  specimen  and  across  its  length  was  found 
to  be  less  than  10F  (~6C) . 

A block  diagram  of  the  instrumentation  used  is 
shown  in  figure  (A4.5).  The  piezoelectric  accelero- 
meter and  charge  amplifier  converted  the  mechanical 
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vibration  to  a voltage  signal  proportional  to  the 
vibration  amplitude.  Electronic  noise  was  filtered 
out  with  the  band-pass  filter  and  the  signal  was 
displayed  on  an  oscilloscope  to  monitor  the  wave 
shape  ^nd  indicate  the  vibration  amplitude.  The 
frequency  counter  indicated  the  vibration  frequency. 

The  level  recorder  gave  the  decay  curve  to  a logarithmic 
scale,  the  slope  of  which  was  used  to  determine  the 
loss  coefficient. 

Before  the  damping  tests  were  performed,  the 
stress  level  on  the  surface  of  a typical  specimen 
was  determined  as  a function  of  the  displacement  at 
the  trigger  so  that  the  stress  amplitude  at  which 
each  test  was  conducted  could  be  monitored  on  the 
oscilloscope.  For  this  calibration,  strain  gages 
were  mounted  on  both  sides  of  the  test  section  of 
the  specimen,  a force  was  applied  to  the  arm  and 
the  displacement  at.  the  trigger  and  the  strain  in 
the  specimen  were  both  recorded  as  shown  in  figure 
(A4.6)  . 

In  order  to  make  measurements  at  170  Hz,  the 
rigid  arm  at  the  free  end  of  the  specimen  was  replaced 
by  a light  rod  as  shown  in  figure  (A4.7).  An 
accelerometer  was  fastened  to  the  end  of  the  rod  to 


measure  the  vibration  amplitude.  A magnetic  coil 
excited  by  an  oscillator  and  audio  amplifier  was  used 
to  drive  the  system  at  its  fundamental  natural 
frequency.  The  frequency  of  the  excitation  was 
varied  until  steady  state  resonant  vibration^ detected 
from  the  oscilloscope  tracef  was  reached.  The  power 
to  the  driving  coil  was  then  cut  off  and  the  system 
allowed  to  decay  freely.  At  lower  frequencies, 
in  reducing  the  data,  it  could  be  assumed  that  all 
of  the  strain  energy  was  in  the  test  specimen.  For 
this  configuration,  the  strain  energy  was  not 
confined  to  the  specimen  alone,  but  some  bending 
strain  energy  was  stored  in  the  relatively  slender 
rod  as  well.  For  a given  amount  of  energy  dissipation 
in  the  specimen  the  measured  loss  coefficient  for  the 
system  will  be  lower.  Computations  based  on  measured 
amplitudes  and  frequencies  showed  that  the  loss 
coefficient  measured  at  this  frequency  had  to  be 
corrected  by  the  factor  shown  ir.  equation  (3.18). 

Natural  Frequency  and  Stress  Analysis 

The  maximum  stress  on  the  surface  of  the  beam 
used  in  the  bending  tests  can  be  calculated  by  the 
method  given  by  Lee  [8].  A sketch  of  the  system 
configuration  is  shown  in  figure  (A4.8).  The  mass 


of  the  beam  itself  is  negligible  compared  with  the 
mass  of  the  mass  at  the  free  end. 


A free  body  diagram  of  the  cantilever  beam  is 
shown  in  figure  (A4.9).  The  deflection  and 
rotation  at  the  end  of  the  cantilever  beam  are 
functions  of  the  bending  moment,  the  shear  force  and 
the  axial  force  at  the  end.  For  small  deformations, 
the  effect  of  the  axial  force  may  be  neglected  and 
the  static  equations  are 


y allV0  + a12M0 
9 = a12V0  + a22M0 

where  the  influence  coefficients  are 

= L^/3EI  = 4.5  * 10  ^ inches/pound 

2 -4 

a^2  = L /2EI  = 4.74  x 10  radians/pound 


(A4.1) 


-4 


a 22  = E/EI  = 6.31  x 10  radians/inch  pound. 


L = 1.5  inches  is  the  length  of  the  beam 

E = 28.5  x 10  psi  is  the  Young's  modulus  of 

Hastelloy  X at  room  temperature 
-4  4 

I = 0.912  x 10  (inches)  is  the  mean  area  moment 
of  inertia  of  the  cross  section  of  the 
specimen. 

The  influence  coefficients  for  the  tapered  beam 
were  measured  for  comparison  with  the  calculated 
values  [8] . The  measurements  were  made  first  by 
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applying  a vertical  force  Vg  at  the  free  end  of  the 
cantilever  beam  (point  A in  figure  (A4.8));  the 
vertical  deflection,  y,  at  point  A and  the  axial 
displacement , x,  at  point  B were  recorded.  The 
coefficients  and  ct^2  were  then  calculated  from 

-4 

ctj-^  = y/Vg  = 6.2  ^ 10  inches/pound 


a12  = (x/*>/vo  = °/V0  " ^ ^ x radians/pound, 

The  coefficients  and  a22  were  found  by  applying 

a horizontal  force  F at  the  point  B and  measuring 
the  vertical  and  horizontal  displacements,  y and  x, 
at  the  same  point.  Then, 

= y/(F£)  = 5.7  x 10  ^ radians/pound 


a22  = (x/i)/(FJ l)  = 9/M0 

-4 

= 6.7  x 10  radians/inch  pound. 

The  measured  coefficients  were  different  from  the 
calculated  ones  probably  because  the  end  conditions 
were  not  ideal.  The  measured  values  were  used  in 
the  calculation  of  the  natural  frequency  and  the 
stress  variation  along  the  surface. 

The  free  body  diagram  for  the  extended  mass 
is  shown  in  figure  (A4.10).  The  equations  of  motion 


for  the  extended  mass  are: 
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M22  = -a12am  - a22  (J  + a m + b m) 


The  natural  frequency  of  the  system  is  found 


from  the  characteristic  equation 
2. 


J-oTM  + K|  = 0 (A4.5) 

where  w is  the  circular  natural  frequency  of  the 
system. 


For  the  test  assembly  at  47  Hz,  we  have 


m = 1.435  x 10  3 lb.  sec2/inch 


-3  2 

J = 8.365  x 10  lb.  sec  inch 


a = 0.85  inch 
b = 1.805  inches 

Using  these  values  and  the  measured  influence 
coefficients,  the  characteristic  equation  of  the 


system  is 


oj4  - (5.5  x 106)  a)2  + (4.72  x 1011)  = 0 


The  solutions  are 


= 2.98  x 10  radians/sec. 


and 


w2  = 2.36  x 10  radians/sec. 


which  gives 


f^  = o)^/  ( 2 7T ) = 47.6  Hz 


and 


f2  = a)2/^2Tr^  = 376  Hz. 


The  measured  value  of  the  lowest  frequency  was  about 
48-49  Hz  for  this  mass. 
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Under  free  vibration,  the  deflection  and  the 
rotation  corresponding  to  the  fundamental 


frequency  are 


iw.t  iu)..t 

y = Y0e  and  6 = eQe 


(A4.6) 


The  fundamental  mode  shape  can  be  found  from 


[-w  M + K]  [XQ]  = 0 


(A4.7) 


where 


[XQ]  = 


*0 


The  solution  of  equation  (A4.7)  gives 


y q/0q  = 0.86  inch 


(A4.8) 


Substituting  (A4.6)  into  (A4.4), 


2 i“lt 

V„  - mu.1  (Y0  + 80)  e 


2 2 2 2 i^i  t 

Mg  = [arnw^  Yq  + (J  + a m + b mJco^O^le 


(A4.9) 


or 


Vq/Mq  = 0.162 


(A4.10) 


The  bending  moment  at  any  cross  section  along  the 
cantilever  beam  (figure  (A4.9))  is 


Mx  = M0  + V0  (L‘X) 


Mx  = MQ  (1  + 0.162  (L-x) 


Substituting  from  equation  (A4.10) 

(A4.ll) 

The  normal  stress  at  any  cross  section  along  the  beam 
is 


ox  - <Mxy)/I 


1 
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where  is  given  by  equation  (A4.ll) 

I = wh3/12 
where 

w = width  of  the  beam 
h = thickness  of  the  beam 
The  stress  on  the  surface  of  the  beam  is  then 
0X  = (6  MQ  (1  + 0.162  (L-x) ) / (wh2) 

= C (1  + 0.162  (L-x) )/h2 
where  C is  a constant. 

At  the  free  end  of  the  beam,  x = 1.4  inches,  and  h = 
0.088  inches.  Therefore, 

Ox  = C(1  + 0.162  x 0.1)/(0.088) 2 
= ±31  C 

At  the  fixed  end,  at  x = 0.1  inches,  h = 0.095  inches, 
Therefore, 

Ox  = c (1  + 0.162  x 1 . 4 ) / ( 0 . 095) 2 
= 136  C 

Thus,  a linear  taper  from  0.095  inches  to  0.088 
inches  over  a gage  length  of  1.3  inches  gives  a 
surface  stress  at  the  fixed  end  which  is  4%  greater 
than  that  at  the  free  end. 
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Appendix  5.  Hollow  Torsion  Cylinder  Apparatus 


This  apparatus  shown  in  figure  (A5.1)  was  used  in 
the  frequency  range  165  Hz  to  600  Hz.  To  obtain 
these  high  frequencies,  the  stiffness  of  the  specimen 
had  to  be  relatively  high.  At  the  same  time,  the 
polar  moment. of  inertia  had  to  be  low  enough  so 
that  measureable  damping  could  be  obtained  with  a 
thin  coating.  The  test  specimen,  shown  in  figure 


(A5.2),  made  of  Hastelloy  X was  in  3 sections,  with 
a total  length  of  15  inches  when  assembled.  The 
hollow  cylindrical  section,  7 inches  long  with  outer 
diameter  0.975  inch  and  wall  thickness  0.050  inch, 

\ < 

formed  the  test  section.  The  middle  4 inches  of 
this  section  was  coated  on  the  outside  with  the 
porcelain  enamel  under  test.  The  ends  of  the 
specimen  fitted  into  the  test  section  by  expansion 
joints  which  could  be  tightened  with  the  tapered 
plug.  Each  end  piece  had  a square  section  for  a 
length  of  1.75  inches  to  permit  clamping.  Only  the 
middle  section  need  be  changed  to  test  various 
coatings.  A torque  wrench  was  used  to  ensure  a 
maximum  torque  of  70  inch  pounds  to  tighten  the  nuts 
which  drew  in  the  tapered  plug,  in  order  to  maintain 
consistent  conditions  for  all  specimens.  The 
specimen  was  horizontal  and  clamped  at  one  end  to  a 
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massive  steel  base.  A light  steel  disc,  with  radial 
arms  diametrically  opposed,  was  clamped  to  the  free 
end  of  the  specimen.  One  of  the  arms  carried  a 
piezoelectric  accelerometer  at  its  end.  The 
excitation  from  a magnetic  coil  was  applied  to  the 
end  of  the  other  arm.  The  excitation  frequency  was 
varied  until  steady  state  resonant  vibration  was 
reached.  The  system  was  allowed  to  decay  freelv  hv 
cutting  off  the  power  to  the  driving  coil  and  the 
decay  rate  was  measured.  The  assembly  shown  in 
figure  (A5.2)  had  a natural  frequency  of  600  Hz.  The 
fundamental  frequency  was  changed  by  clamping  steel 
weights  to  the  free  end  of  the  specimen.  The 
frequencies  so  obtained  were  430  Hz,  340  Hz  and  165  Hz. 

A cylindrical  electric  furnace  was  used  to  heat 
the  specimen.  The  furnace  had  a bore  of  2.5  inches 
and  was  10  inches  long  so  that  2.5  inches  of  the 
specimen  projected  outside  the  furnace  at  either 
end  for  clamping,  measurement  and  control.  Insulating 
plates  were  used  at  either  end  to  close  the  furnace 
chamber.  The  furnace  winding  was  tapped  at  a number 
of  points  and  brought  out  to  terminals  on  the  outside 
of  the  furnace.  Certain  sections  of  the  furnace 
winding  could  be  shunted  to  maintain  a uniform 
temperature  distribution  in  the  furnace.  Three  chromel- 
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alumel  thermocouples,  their  lead  wires  enclosed  in 
refractory  tubes,  could  be  inserted  into  the  furnace 
chamber  through  holes  in  the  wall  of  the  furnace. 

They  covered  a length  of  about  3 inches  in  the  middle 
of  the  furnace  and  could  be  radially  adjusted  to  be 
about  0.03  inch  away  from  the  specimen  surface. 

The  middle  thermocouple  was  used  to  control  the 
furnace  temperature  and  was  also  the  indicating 
thermocouple  during  the  damping  tests.  A temperature 
survey  was  conducted  with  thermocouples  mounted  at 
various  positions  on  the  apparatus  with  an  uncoated 
specimen  as  shown  in  figure  (A5.3).  Thermocouple 
#7  was  located  at  the  position  of  the  accelerometer 
and  a blower  was  used  to  keep  its  temperature  within 
the  safe  operating  range  of  the  accelerometer.  The 
resistance  shunts  across  the  various  sections  of  the 
furnace  winding  were  adjusted,  until  the  temperature 
variation  in  the  middle  4 inches  of  the  specimen 
was  less  than  20F  (11C) . The  indicating  thermocouple 
(#8)  was  found  to  show  a temperature  20F  (11C)  higher 
than  that  of  the  specimen.  A soak  time  test  showed 
that  it  took  at  least  two  hours  to  reach  temperature 
equilibrium  for  the  first,  relatively  large, 
temperature  change,  and  about  45  minutes  for  each 
subsequent  smaller  temperature  change. 
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The  accelerometer  and  a charge  amplifier  were 
used  to  produce  a voltage  signal  proportional  to  the 


i 


vibration  amplitude.  The  signal  was  displayed  on 
an  oscilloscope  to  monitor  the  wave  shape  and  indicate 
the  vibration  amplitude.  The  frequency  \*7as  measured 
with  a counter.  A level  recorder  was  used  to  measure 
the  decay  rate.  The  magnetic  coil  used  to  excite 
the  system  was  driven  by  a low  freqeuncy  oscillator 
and  audio  amplifier. 

A calibration  of  the  stress  level  versus  the 
angular  displacement  at  the  free  end  was  made  on  an 
uncoated  specimen.  45°  strain  gages  for  shear 
measurements  were  mounted  at  either  end  of  the  test 
section  of  the  specimen.  Torque  was  applied  by 
hanging  weights  from  the  end  of  a lever  arm 
fastened  to  the  free  end  of  the  specimen.  A dial 
gage  showed  the  linear  displacement  at  the  end  of  a 
light  aluminum  arm;  for  small  displacements,  this 
reading  divided  by  the  distance  to  the  center  of 
the  specimen  gives  the  angular  displacement.  Figure 
(A5.4)  is  a graph  of  the  average  strain  in  the 
.pecrmen  versus  the  angular  displacement  at  the  free 
end . 

The  lowest  frequency  of  165  Hz  was  obtained  with 
this  apparatus  by  clamping  a solid  steel  disc  shown 
in  figure  (A5.5)  to  the  end  of  the  specimen.  The 
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ioss  coefficient  at  this  frequency  was  found  to  be 
amplitude  dependent.  The  apparatus  was  set  up  so 
that  the  specimen  was  vertical,  to  eliminate  static 
deflection  of  the  specimen  due  to  the  heavy  weiyht 
at  the  end;  the  non-linearity,  however,  still 
remained.  Damping  measurements  at  room  temperature 
were  made  on  a solid  steel  specimen,  0.75  inches 
in  diameter  and  15  inches  long.  This  showed  no 
dependence  of  the  loss  coefficient  on  vibration 
amplitude.  Apparently,  the  expansion  fits  of  the 
test  section  with  the  end  grips  caused  this  non- 
linearity. The  temperature  dependence  of  the  loss 
coefficient  of  the  system  with  a coated  and  an 
uncoated  specimen  was  determined  at  various  stress 
levels.  The  temperature  dependence  of  the  loss 
coefficient  due  to  the  coating  alone  was  determined 
at  various  stress  levels  as  the  difference  in  the 
damping  value  between  the  coated  and  uncoated 
specimens  al  the  same  stress  amplitude.  Figure  (A5-6) 
shows  the  measured  system  damping  for  the  coated  and 
bare  specimens  at  two  different  temperatures  as 
functions  of  amplitude.  Since  the  difference  in  ns 
between  the  coated  and  uncoated  specimens  is 
independent  of  amplitude,  the  damping  in  the  coating 
alone  is  not  amplitude  dependent;  the  non-linearity 
is  due  entiiely  to  mechanical  misfit  in  the  expansion 


joints.  Al]  damping  values  for  this  frequency  were 
determined  at  an  amplitude  corresponding  to  about 
0.1  volts  on  the  abscissa  in  figure  (A5.6)  which 
corresponds  to  a shear  stress  level  of  50  psi  on 
the  surface  of  the  specimen.  Similar  amplitude 
surveys  were  made  with  other  discs  at  higher 
frequencies  and  no  such  non  linearity  was  observed. 


Natural  Frequency 


The  test  specimen  is  shown  in  Figure  (A5.2). 

It  was  held  at  one  end  by  a clamp  1.75  inches  long, 
so  that  the  effective  length  of  the  specimen  was 
13.25  inches.  The  stiffness  distribution  in  the 
specimen  was  used  to  calculate  the  strain  energy 
for  data  reduction.  The  measured  and  computed 
natural  frequencies  were  compared  to  verify  the 
strain  distribution  in  the  test  section. 

The  Holzer  method  [47]  was  used  for  the  frequency 
calculation.  Experience  shows  that  calculated  mass 
polar  moments  of  inertia  are  quite  accurate,  but 
that  torsional  stiffnesses  are  more  difficult.  The 
stiffnesses  of  the  specimen,  between  the  support  and 
three  different  points  on  the  specimen,  were 
measured  by  applying  a known  torque,  and  measuring 
the  corresponding  angular  displacement.  These 
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points,  A,  B and  C were  at  4.25  inches,  7.25  inches 

and  12.50  inches  respectively  from  the  support. 

Figure  (A5.7)  shows  the  graph  of  the  measured  angular 

displacement  at  these  points  versus  the  applied 

torque.  The  slopes  in  this  graph  give  the  stiffnesses 

between  the  support  and  the  corresponding  points. 

The  stiffnesses  K, , K„  and  K„  were  also 

A B C 

calculated  from  the  polar  moment  of  inertia 

Now, 

1/KA  = VM  * 5 Vip  a* 

where  G is  the  shear  modulus,  and  I is  the  polar 
moment  of  inertia.  Then: 

1/K.  = 1/G  97.04  = 1/(11  X 106)  X 97.04 

/i 

= 8.73  x 10  ® radians/inch  pound 
or 

4 

KA  = 11.45  x 10  inch  pound/radian 

4 

Similarly,  Kg  = 5.95  x 10  inch  pound/radian , Kc  = 

3.56  x 104  inch  pound/radian.  The  measured  stiffnesses, 

with  superscript  m,  from  Figure  (A5.7)  are 

m 4 m 4 

Ka  = 6.0  x 10  inch  pound/radian,  Kg  =3.96x10 

inch  pound/radian  and  K^,1'1  = 2.13  x 104  inch  pound/ 

radian.  The  measured  stiffnesses  in  each  case  were 

about  0.6  times  the  stiffness  calculated  from  the 

polar  moment  of  inertia.  The  difference  was  probably 

due  to  the  end  support  not  being  perfectly  rigid  and 
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some  loss  of  stiffness  in  the  expansion  joints  between 
the  test  section  and  the  end  grips.  In  the  calculation 
of  the  natural  frequency,  the  stiffnesses  obtained 
from  the  polar  moment  of  inertia  were  used  after 
correcting  them  by  the  factor  of  0.6.  The  specimen 
was  divided  into  53  segments  each  of  length  0.25 
inches.  The  mass  moment  of  inertia  of  the  driving 
arm  was  added  tc  that  of  disc  #1.  Figure  (A5.8)  shows 
the  graph  of  assumed  frequency  versus  the  residual 
angular  displacement  at  the  clamped  end.  The  natural 
frequency  of  the  system  can  be  read  off  as  644  Hz. 
Figure  (A5.9)  shows  the  mode  shape  at  that  frequency. 
The  measured  frequency  of  the  specimen  was  605  Hz. 

The  lower  frequency  of  170  Hz  was  obtained  by 
fastening  a heavy  disc  near  the  free  end  of  the 
specimen,  at  the  point  C.  Dividing  the  specimen 
into  N parts  as  before,  the  total  kinetic  energy  of 

the  system  is 

2 N 2 2 
KE  = a>  ( Z J.  0 + J 0 )/2 

i=1  ii  mm" 

where  J = mass  moment  of  inertia  of  the  heavy  disc 
m 

and  = mass  moment  of  inertia  of  the  discs  of  the 
Holzer  division.  0's  are  the  corresponding  angular 
displacements  and  w = frequency  in  radians/sec. 
includes  the  driving  arm;  but  each  of  the  J^'s 
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for  i 


1 


- 2 .... N is  negligibly  small  compared  to  J 


m 


Therefore , 


KE  = co2  (J  62  + J 02)/2 
11  mm 


Then  for 


0m  = 1,  0-j^  = (13 . 25)  /12 . 5 = 1.06 

— 3 2 

= 1.0  x io  inch  pound  second 

- 3 2 

Jm  = 19  x 10  inch  pound  second 
Therefore , 

2 

KE  = uj  x 10.1  x io  inch  pound 
The  potential  energy  = Kc0m/2. 

4 

K^,  measured  = 2.13  x io  inch  pound/radian 
Therefore , 


PE  = 1.07  x io4  inch  pound 
Equating  kinetic  and  potential  energies, 
w = 1.06  x io®  (radians)2 
Therefore , 

frequency  f = 166  Hz. 

The  measured  frequency  with  this  mass  was  175  Hz. 
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Appendix  6.  Details  of  Experimental  Procedures 


This  appendix  describes  the  various  procedural 
problems  that  were  encountered  during  the  experiments 
and  the  measures  taken  to  overcome  them. 

In  order  to  reduce  variability  in  system 
measurements  from  test  to  test,  consistent  mounting 
conditions  had  to  be  maintained  for  all  specimens 
and  extraneous  damping  had  to  be  reduced  to  a 
minimum.  A torque  wrench  was  used  to  tighten  bolts 
and  nuts  to  the  same  torque.  In  the  bending 
apparatus,  the  high  temperature  bolts  used  were 
tightened  to  a torque  of  120  inch  pounds,  and  those 
in  the  high  frequency  torsion  apparatus  to  70  inch 
pounds.  In  assembling  the  latter  specimen  a jig, 
exactly  15  inches  long  was  used  so  that  the  end 
grips  could  be  positioned  correctly  inside  the 
center  hollow  cylinder.  In  both  the  torsion 
experiments  the  weights  and  the  driving  arm  had  to 
be  fastened  at  the  free  end  of  the  specimen  at  the 
same  point  from  test  to  test.  In  the  high  frequency 
torsion  experiment  the  driving  arm  was  placed  at 
the  very  end  of  the  specimen.  The  clearance  between 
the  insulating  plates  enclosing  the  furnace  and  the 
driving  coil  support  was  only  an  inch;  the  weights 
were  positioned  in  the  middle  of  this  clearance. 
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In  the  torsion  pendulum  apparatus,  the  position  of 
the  weight  and  driving  arm  was  determined  by  a 
pair  of  cylinders  placed  on  the  fixed  floor  of  the 
apparatus . 

To  reduce  the  room  temperature  damping  of  the 
systems,  they  were  isolated  from  the  building 
vibration  using  rubber  pads  (isomode  pads) . The 
various  parts  of  the  apparatus  had  to  be  properly 
aligned.  The  ball  bearing  in  the  torsion  pendulum 
case  was  particularly  difficult  to  align.  The 
slightest  misalignment  there  caused  the  system 
damping  at  room  temperature  to  increase  from  0.0002 
to  0.001.  The  cable  from  the  accelerometer  was 
taped  to  the  base  after  adjusting  the  tension  and 
position  of  the  cable  to  give  minimum  damping. 

Room  temperature  measurements  were  always  made  with 
the  furnace  chamber  closed  with  the  insulating  plates 
so  that  variations  in  air  friction  between  room 
temperature  and  higher  temperatures  were  minimized. 
The  lowest  room  temperature  damping  measured  was 
about  0.00015;  the  highest  was  over  0.002  for  the 
high  frequency  torsion  experiment.  The  high 
frequency  electronic  noise  from  the  instrumentation 
used  and  the  line  frequency  noise  were  filtered  out 
using  a band-pass  filter. 
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All  temperature  measurements  were  made  with 


chrome 1 -a lumel  thermocouples.  The  same  thermocouple 
was  used  to  indicate  the  specimen  temperature  for 
about  10  tests,  after  which  the  thermocouple 
junction  was  renewed.  The  drift  in  the  indicated 
temperature  due  to  the  on-off  furnace  cycle  was 
about  1C.  In  the  torsion  experiments  the  indicating 
thermocouple  was  not  in  contact  with  the  specimen, 
but  about  0.03  inch  away  from  the  surface  of  the 
specimen.  Before  the  damping  tests  were  conducted, 
a temperature  survey  was  done  using  several  thermo- 
couples as  shown  in  figures  (A3. 4)  and  (AS. 3), 
fastened  along  the  length  of  specimen.  These 
thermocouples  had  their  junctions  in  good  contact 
with  the  specimen.  The  temperatures  indicated  by 
these  thermocouples  were  compared  with  the  temper- 
ature indicated  by  the  thermocouple  used  in  the 
damping  tests  at  a number  of  different  temperatures 
in  the  range  of  interest.  It  was  found  that  the 
temperature  indicated  during  the  tests  was  liC 
lower  for  the  torsion  pendulum  apparatus  and  11C 
higher  for  the  high  frequency  torsion  apparatus  when 
the  furnace  resistances  were  set  to  give  uniform 
temperature  across  the  4 inch  test  section  of  the 
specimens.  During  the  tests,  sufficient  time  was 
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allowed  for  the  specimen  to  attain  temperature 
equilibrium  - 2 hours  for  the  first  temperature 
change  and  4 5 miniates  for  each  subsequent  temperature 
change . 

The  accelerometer  used  in  the  earlier  tests 
had  a maximum  safe  operating  temperature  of  180C. 

It  had  to  be  kept  within  that  limit  by  blowing  air 
over  it.  The  blower  was  used  even  during  the 
temperature  surveys.  It  lowered  the  temperature 
of  the  specimen  by  about  3C.  For  some  later 
experiments,  a high  temperature  accelerometer 
capable  of  operating  safely  up  to  4C0C  was  used 
and  the  blower  was  then  not  required. 

In  order  to  avoid  geometry  changes  due  to 
differences  in  thermal  expansion,  all  parts  of  any 
apparatus  subjected  to  high  temperatures  were  made 
of  the  same  material.  High  temperature  bolts  and 
nuts  were  used.  To  prevent  binding  of  screw  threads 
exposed  to  high  temperatures,  a coating  of  ceric 
oxide  was  used  on  them.  Care  was  taken  not  to 
exceed  a torque  of  70  inch  pounds  in  tightesing 
the  high  temperature  bolts  which  draw  in  the  tapered 
plug  in  the  high  frequency  torsion  experiment.  A 
higher  torque  at  room  temperature  did  not  change 
the  system  damping  appreciably;  too  high  a torque 
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might  cause  the  different  parts  of  the  specimen  to 
weld  together  when  heated  to  temperatures  exceeding 
1000C. 


The  frequency  measurements  were  made  with  a 
counter  whose  gate  could  be  set  for  a 1,  10  or  100 
second  measurement.  For  very  low  frequencies  as  in 
the  torsion  pendulum  experiments,  the  decay  was 
sufficiently  slow  that  the  measurement  could  be 
made  over  a 100  second  interval.  However,  the 
chances  of  stray  noise  signals  being  picked  up  and 
counted  were  also  increased.  The  count  on  the 
instrument  had  to  be  watched  carefully  over  the 
entire  measurement  interval  to  ensure  that  no  stray 
signals  were  counted.  By  taking  an  average  of 
several  such  frequency  measurements,  the  error  in 
the  frequency  at  1.2  Hz  could  be  controlled  to 
0.2%.  For  the  frequencies  170  Hz  and  above,  since 
a magnetic  driver  was  used  to  excite  the  system, 
the  resonant  frequency  could  be  measured  before 
cutting  off  the  power  to  the  driver.  By  making  a 
10  second  measurement,  the  error  was  limited  to  a 
maximum  of  0.1%.  When  using  a magnetic  driver,  the 
output  of  the  oscillator  and  the  signal  from  the 
accelerometer  were  fed  to  the  horizontal  and  vertical 
terminals  of  the  oscilloscope  respectively.  The 


resonant  frequency  of  the  system  was  determined  from 
the  Lissajous  figure  displayed  on  the  oscilloscope. 

To  ensure  that  the  system  was  being  driven  exactly 
at  resonance  before  making  the  freqeuncy  measurement, 
the  system  was  allowed  to  decay  freely  by  cutting  off 
the  power  to  the  oscillator.  The  collapse  of  the 
Lissajous  figure  on  the  oscilloscope  was  watched  for 
any  phase  changes  indicated  by  a change  in  shape  of 
the  figure.  The  frequency  of  the  driver  was 
adjusted  so  that  the  Lissajous  figure  decayed 
vertically  without  changing  shape.  The  least 
accurate  of  the  frequency  measurements  were  the  ones 
made  at  47  Hz.  The  decay  rate  in  the  transition 
region  was  too  high  to  permit  a 10  second  measurement. 
By  taking  an  average  of  about  30  measurements  over 
a 1 second  interval,  the  error  could  be  controlled 
to  0.3%. 

The  thickness  of  the  coating  could  not  be 
controlled  accurately  during  application.  It  had 
therefore  to  be  measured  after  application.  It  was 
obtained  as  the  difference  in  specimen  thicknesses 
(diameters)  before  and  after  application  of  the 

-4 

coating.  A micrometer  with  a least  count  of  10 
inch  was  used  for  the  torsion  specimens.  However, 
for  the  bending  specimens,  a micrometer  was  not  found 
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satisfactory  because  of  the  taper.  A jig  was  devised 
to  measure  the  thickness  of  the  bending  specimen. 

It  consisted  of  a plate  with  V-groves  exactly  1.3 
inches  apart  which  was  the  gage  length  of  the  bending 
specimen.  Two  roller  pins,  with  radius  slightly 
less  than  the  radius  of  the  fillets  at  the  two  ends 
of  the  test  section  of  the  specimen,  were  placed  one 
in  each  groove.  The  specimen  was  then  laid  over  the 
i pins  so  that  the  lower  surface  of  the  test  section 

was  horizontal  at  the  level  of  the  top  of  the  pins. 

A dial  gage  with  a least  count  of  10  4 inch  was 
used  to  measure  the  specimen  thickness,  the  reading 
at  the  top  of  the  two  pins  was  the  same  and  set  to 
zero.  The  specimen  thickness  was  then  read  off  the 
dial  gage  at  various  points  on  the  specimen.  Care 
was  taken  to  see  that  the  readings  were  taken  at  the 
same  points  with  both  coated  and  uncoated  specimens. 
For  the  torsion  specimens  also,  the  thickness  was 
measured  very  precisely  at  several  points  and  an 
average  taken.  The  error  in  the  specimen  thickness 
so  obtained  is  not  expected  to  exceed  0.05%. 


TABLE  (2.1) 


COMPOSITION  OF  ENAMEL  FRITS 


T 

C.V. #12325,% 

C.V. #16845,% 

Silicon  Dioxide 

(Si02) 

l 

33.68 

27.46 

Aluminum  Oxide 

(ai2o3) 

3.25 

3.09 

Iron  Oxide 

(Fe2°3} 

0.24 

0.28 

Sodium  Oxide 

(Nao0)  1 

13.27 

13.63 

Calcium  Oxide 

(CaO) 

5.19 

5.47 

Zinc  Oxide 

(ZnO) 

3.11 

5.69 

Barium  Oxide 

(BaO) 

12.72 

12.32 

Potassium  Oxide 

(K20) 

18.64 

18.85 

Boron  Oxide 

(B2°3) 

8.50 

8.50 

Magnesium  Oxide 

(MgO) 

0.15 

1.01 

Lead  Oxide 

(PbO) 

0.02 

i 

0.02 

J 

Antimony  Oxide 

(SbO) 

1 

1.76 

i 

Manganese  Oxide 

(Mn02) 

i 

0.01  j 

Molybdenum  Oxide 

(Mo02) 

- 

0.74 

Trace  Constituents 

Antimony 

Titanium 

Copper 

Manganese 

Silver 

Titanium 

! 

Chromium 

1 

..... 
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TABLE  2.3 


VALUES  OF  2f (x ' , 8)  = 2M" (x’)/ 


TABLE  4.1 

MANUFACTURERS  OF  COATING  FRITS 


Type  of  Coating 

Firing  Temperature,  °C 

Manufacturer 

AL-2 (A) 

550 

Ferro 

S5-8A 

550 

Solar 

CV-145-S (B) 
or  #16845 

600 

Chicago 

325 (Bb) 
or  #12325 

730 

Chicago 

S9919I 

900 

Solar 

S6100M 

980 

Solar 

Addresses : 
Ferro : 

Solar : 

Chicago : 


Ferro  Corporation 
One  Erieview  Plaza 
Cleveland,  Ohio  44114 

Solar,  Division  of  International  Harvester  Company 

2200  Pacific  Highway 

San  Diego,  California  92112 

Chicago  Vitreous  Corporation 
1425  South  55th  Court 
Cicero,  Illinois  60650 


MAXIMUM  DAMPING  AND  TEMPERATURE  DATA  - TORSION  PENDULUM 


0.0011  | 682  0.0088 


TABLE  4.4 


MAXIMUM  DAMPING  AND  TEMPERATURE  DATA  - 
HIGH  FREQUENCY  TORSION 

SPECIMEN  #1 

COATING  TYPE:  #12325 

COATING  THICKNESS:  0.0057  INS 


165  HZ 

340  HZ 

430  HZ  i 600  HZ 

1 

T,°C  n 
* c max 

j 

T°c  n 

’ c max 

• 

T,°C  n ' T,°C  n 

5 c max  * c max 

587  0.0054 

599  0.0056 

i 

! | 

,599  0.0056  ,601  0.003  1 

|i  ! i 
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TABLE  4.5 


VALUES  OF  G2'  AND  G2"  - TORSION  PENDULUM 


TABLE  4 . 6 


VALUES  OF  E^  AND  E^  - BENDING 


— 

SP.# 

Coating 

Type 

— 
Coating 
Thickness 
d2,  ins.  ; 

E'  psi 

* ID'6 
at  Room  T. 

E2pfj ; 

X 10  6 
Maximum 

5 

S9919D 

0.0024 

10.0 

1.75 

7 

ft 

0.0057 

9.5 

1,84 

16 

#12325 

0.00185 

13.5 

1.98 

1 

16A 

II 

0.0029 

10.0 

1.34 

1 

22 

11 

0.004 

8.2 

t— 1 
* 

O 

00 

l5 

#16845 

0.002 

9.0 

1 

1.28 

15A 

II 

0.0023 

9.0 

1.30 

Note:  The  modulus  values  have  been  corrected  for 

anticlastic  curvature  by  multiplying  the 
measured  values  by  0.9.  This  corresponds  to 
a Poisson's  ratio,  v,  of  0.31. 
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Coated  at  the  University  of  Minnesota 


TABLE  4.9 
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TEMPERATURE 

FIG.  2.1.  VOLUME  VS  TEMPERATURE  FOR  GLASS . 


TEMPERATURE  C 


FIG.  2.2.  E2'  AND  E2"  VS  TEMPERATURE 

COATING  S9919D,  SP  #5,B. 


VS  TEMPERATURE;  f = 47  HZ,  SP 
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FIG.  4.9.  n£  VS  TEMPERATURE  FOR  DIFFERENT 
FREQUENCIES;  SP.  #2,  TP 


TEMPERATURE  C 


FIG.  4.11.  r\  VS  TEMPERATURE  FOR  DIFFERENT 
c 

FREQUENCIES;  SP.  #1,  T. 
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FIG.  4.17.  E2"  VS  TEMPERATURE  FOR  DIFFERENT 

FREQUENCIES;  SP.  #16A,  B. 
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FIG.  4.19.  G2"  VS  TEMPERATURE;  COATING  #12325 
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APPARATUS  FOR  TORSION  PENDULUM  TESTS 


'DECKER  corporation  band-pass  COUNTER - 

1906  DELTA  UNIT  1 FILTER  TIMER 


(1)  FURNACE 

(2)  SPECIMEN 

(3)  MICROMETER 

(4)  ACCELEROMETER 

(5)  STEEL  PLATE 


(1)  SPECIMEN 
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